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ABSTRACT 

We have analyzed photometry from space- and ground-based cameras to iden- 
tify a// bright red giant branch (RGB), horizontal branch (HE), and asymptotic 
giant branch (AGB) stars within 10' of the center of the globular cluster M13. 

We identify a modest (7%) population of HB stars redder than the primary 
peak (including RR Lyrae variables at the blue end of the instability strip) that is 
somewhat more concentrated to the cluster core than the rest of the evolved stars. 
We find support for the idea that they are noticeably evolved and in the late stages 
of depleting helium in their cores. This resolves a disagreement between distance 
moduli derived from the tip of the red giant branch and from stars in or near 
the RR Lyrae instability strip. We identified disagreements between HB model 
sets on whether stars with Tefj < 10000 K (near the "knee" of the horizontal 
branch in optical CMDs) should evolve redward or blueward, and the differences 
may depend on the inclusion of diffusion in the stellar interior. The sharp cut at 
the red end of M13's HB provides strong evidence that stars from the dominant 
HB group must still be undergoing blue loops, which implies that diffusion is 
being inhibited. We argue that M13's HB is a somewhat pathological case — 
the dominant HB population occurs very near the "knee" in optical CMDs, and 
evolved stars exclusively appear redward of that peak, leading to the incorrect 
appearance of a continuation of the unevolved HB. 

We identify two stars as "blue hook" star candidates — the faintest stars 
in optical bands that remain significantly subluminous in the shortest ultravi- 
olet wavelength photometry available. M13 also has a distinct group of stars 
previously identified with the "second U jump". Based on far UV photome- 
try, we find that these stars have genuinely high temperatures (probably 26000 
K ^ < 31000 K), and are not produced by a jump in brightness at lower 
temperature (Tgfj ~ 22000 K) as previously suggested. These stars are brighter 



- 2 - 



than other stars of similar color (either redder or bluer), and may be examples 
of "early hot flashers" that ignite core helium fusion shortly after leaving the red 
giant branch. We used ultraviolet photometry to identify hot post-HB stars, and 
based on their numbers (relative to canonical AGB stars) we estimate the posi- 
tion on the HE where the morphology of the post-HB tracks change to / ~ 17.3, 
between the two peaks in the HE distribution. 

Concerning the possibility of helium enrichment in M13, we revisited the 
helium-sensitive R ratio, applying a new method for correcting star counts for 
larger hfetimes of hot horizontal branch stars. We find that MlS's R ratio is in 
agreement with theoretical values for primordial helium abundance Yp — 0.245 
and inconsistent with a helium enhancement Ay = 0.04. The brightness of 
the horizontal branch (both in comparison to the end of the canonical HE and 
to the tip of the red giant branch) also appears to rule out the idea that the 
envelopes of the reddest HE stars have been significantly enriched in helium. The 
absolute colors of the turnoffs of M3 and M13 may potentially be used to look 
for differences in their mean helium abundances, but there are inconsistencies 
in current datasets between colors using different filters that prevent a solid 
conclusion. 

The numbers of stars on the lower red giant branch and in the red giant 
bump agree very well with recent theoretical models, although there are slight 
indications of a deficit of red giant stars above the bump. There is not convincing 
evidence that a large fraction of stars leave the RGE before undergoing a core 
helium fiash. 

Subject headings: stars: AGB and post- AGB — stars: evolution — Hertzsprung- 
Russell and C-M diagrams — stars: horizontal branch — globular clusters: in- 
dividual (M13) — ultraviolet: stars 



1. Introduction 

M13 was one of the first globular clusters identified as having unusually blue horizontal 
branch (HB) stars, and it remains one of the prototypes of the "long blue tail" with stars 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Science Institute. STScI is operated by the Association of Universities for 
Research in Astronomy, Inc. under NASA contract NAS5-26555. 
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approaching the main sequence for hehum stars. As ISmithI ( 120051 ) notes, differences between 



Barnard 


( 


1909. 


1914) 



httle reddened cluster M3, M13 forms half of the best known "second parameter" pair. 
The chemical composition of the hydrogen envelopes of HB stars affects their observable 
properties via opacity and mean molecular weight. The "first parameter" is heavy element 
content, where higher metallicity produces higher envelope opacity and generally redder 
stars. Because M3 and M13 have near ly the same iron ab undances (([Fe/H]) = —1.53 for 
M13 versus ([Fe/H]) = -1.45 for M3; ISneden et all \20oi \ but M3 has a much redder HB 



including a huge number of RR Lyrae variable stars, a second parameter is needed. In 
addition to the color shift between the HBs of M3 and M13, the HB stars in M13 show a 
bimodal distribution that is not present in M3. On its own, this fact implies that the HB 
stars were produced by at least two different populations of cluster stars or involve different 
methods of producing HB stars. M3 and M13 share a number of similarities beyond iron 
abundance, and we tabulate some of their characteristics in Table [H 

Many theories have been proposed to explain the HB differences, and two of the main 
goals of this paper are to 1) assemble a large and complete set of photometric data for 
M13, and 2) use the photometric data to examine questions bearing on the production of 
M13's horizontal branch stars. Because of the complexity of the HB, it is very doubtful 
that one explanation can cover all of its aspects. Before we describe our results, we briefly 
summarize the main hypotheses we will examining, and the primary reasons they are viable. 
We emphasize that they are not mutually exclusive. 

The At Hypothesis. Early models showed that as the mass of the hydrogen-rich en- 
velope of an HB star is decreased, the surface temperature increases with relatively little 
change in luminosity. In this hypothesis, age differences between populations of stars lead to 
differences in mass between stars leaving the main sequence and between stars reaching the 
HB. While it is natural to expect that clusters in the Milky Way were born at different times, 
age differences a re hard to prove except for clusters that are much younger than the average. 
Rey et al.l ( 1200 ll ) compared M13 with M3 and found that turnoff-to-giant branch color dif- 
ferences and changes i n HB morpholog y were consistent with an age difference At = 1.7 ±0.7 
Gyr (with M13 older). ICatelaru ( 120091 ) finds that the age differences implied by differences in 
the cluster CMDs near the turnoff can explain the HB morphology as long as M3 is younger 
than about 12 Gyr and the EHB stars in M13 are presumed to arise from a process that 
is unrelated to the age. However, there are aspects of M13's population that cannot be 
explained in this hypothesis. For example, neither of the studies above could reproduce the 
bluest HB stars in M13 in synthetic HB simulations using the same chemical composition 
and dispersion in stellar mass used for M3. 
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The Ay Hypothesis. Variations in helium content [Y) result in differences in position on 
the HB, largely because great er helium abundance allow s lower mass stars to leav e the main 
sequence at the present day (ID'Antona fc Caloill2008bl ). Ijohnson fc Boltd ( 1l998l ) proposed 
that a helium abundance difference AY ~ 0.05 (with M13 the more helium rich) could 
be responsible for many of the unusual features of the color-magnitud e diagram (CMP), 
including an interesting difference in the slopes of the subgiant branch. ICaloi fc D'Antona 
(I2OO5I ) also examined data for M3 and M13, finding the luminosity of the red giant bump and 
RR Lyrae stars relative to the cluster turnoff are consistent with an enhancement AY ~ 0.04. 
This picture has been taken very seriously with the discovery of multiple stellar populations 
in some clu sters. For example , NGC 2808 was found to have at least thre e identifiable main 
sequences ( iPiotto et al.l 120071 ). while a; Cen has a blue main sequence (IBedin et al.l |2004[ ) 
that appears to be helium enriched (jPiotto et al.l 120051 ). u Cen and NGC 2808 are among 
the most massive clusters known in the Milky Way, which may enable them to retain gas 
that has been processed and released by a first generation of stars. 

While the helium abundance is very difficult to measure except in limited circumstances, 
spectroscopic observations of other heavy element species lead to the belief that helium 
was probably enriched in some clusters. Stars in M13 are well-known to have star-to-star 
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enrichment can only be acc omplished in hydrogen- fusion regions where significant production 
of helium is accomplished ( iDenisenkov &: Denisenkovalll989l : iLanger et al.lll993l ). and star- 
to-star variations on the main sequence require that they must have been present in the gas 
forming the stars. 

The Ay hypothesis is attractive for M13 because it may explain the blueward shift of 
the main body of HB stars compared to M3's po pulation, and the bimodal ity of the HB 
(as an additional population within the cluster). iD'Antona &: Caloil (l2008b( ) conducted a 
fit to M13's HB using helium-enriched models, and found that a fit required 70% of the 
population to be enriched to 0.27 <Y< 0.35 (a mean AY pa 0.04), and the remaining 30% 
to be enriched to Y ^ 0.38. There are some difficulties with this picture though. In the 



D'Antona fc Caloil models of M13, they still needed to assume a rather large (but constant) 
total mass loss on the RGB (O.I8M0). Because M13 has few HB stars in the instability strip 
or redward (where M3's HB is heavily populated), virtually a//of M13's HB stars must also be 
more helium rich than M3's. This is in striking contrast to more massive clusters that show 
strongly bimodal HB star distributions in which the redder HB population is interpreted 
as a first generation of stars formed from primordial material, while subsequent generations 
have varying degrees of enrichment and are bluer. The massive clusters that are inferred 
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to have such large spreads in Y generally also have multiple main sequence or subgiant 
branch populations, whereas M13 has shown no sign of multiple populations to date. There 
is also not a clear bimodality in the spectroscopic abundances and a definite connection has 



never been made between the abundances and HB morphology (although see ICarretta et al. 



20071 . for indications that the maxir num T^ff ex t ent o f the HB is correlat ed with the exten t 
of observed Na-0 anticorrelations) . iRey et al.l (120011 ) attribute some of I Johnson fc Boltd 's 
conclusions to slight missteps in the implementation of their relative age comparison. 

The M Hypothesis. From early models, it was recognized that a significant amount of 
mass loss is needed, probably on the red giant branch (RGB), to produce the colors of the 
majority HB populations in most clusters. Dispersion in colors was then taken to mean 
that there are star-to-star differences in mass loss, although a mechanism to produce these 
differences has not been identified. 

Independent of the majority of HB stars, there is a population that seems to require 



strong mass loss: the "blue hook" stars ( 


Castellani & Castellani 


1993: D'Cruz et al. 1996: 
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2008 


). In ultraviolet color-magnitude dia- 


grams of some of the most massive clusters {'. 
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are found fainter and redder than the zero-age HB (ZAHB) at its blue end, meaning that 
they must have almost no hydrogen envelope. If a star loses virtually all of its hydrogen 
envelope before reaching the tip of the RGB, it can leave the RGB without igniting core 
helium fusion. As the star contracts onto the He white dwarf cooling curve, a late He fiash 



can b e ignited that drives a convection zone that reaches hydrogen rich layers (IBrown et al. 
200lh . 



The Evolution Hypothesis. As relatively cool HB stars convert He into carbon and oxy- 
gen, they are expected to eventually evolve brightward and redward toward the asymptotic 
giant branch (AGB). Depending on the distribution of stars on the HB, evolving HB stars 
could be mistaken for fainter, more slowly evolving HB stars, thereby misrepresenting the 
brightness of the HB. Clusters with large blue HB populations (like M13) are most suscepti- 
ble to this effect because evolutionary tracks may nearly parallel the ZAHB at the blue end 
of the HB distribution where the relative number of unevolved stars drops rapidly. Because 
the HB is a frequently used standard candle in astronomy, it is worth studying the degree 
to which this affects stellar populations. 

Evolutio nary effects ha v e been inferred from the pulsation properties of RR Lyrae stars. 
For example, lJurcsik et al.l (120031) used magnitudes and peri ods to identify RR Lyraes in 
different stages of their HB evolution. ICacciari et al.l (l2005[ l identified mean lines in the 
period-amplitude diagram for different subsets of M3 variables, and labeled them as regular 
or "well-evolved" . They showed that in at least some other clusters, the majority of variables 
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could be identified with one group or the other. 

For the purposes of this paper we focused on post main- sequence evolution. Using 
datasets from telescopes and instruments having a wide range in spatial resolution and 
field size, we attempted to completely survey the evolved stars from the center far into 
the outskirts of the cluster. We discuss the observational material and the analysis of the 
photometry in §2. In §3, we describe the steps used to identify the evolutionary status of the 
evolved cluster stars. We examine the red giant branch and horizontal branch populations 
in greater detail in §4 and §5, respectively. In ^ we look at population ratios and their 
relationship to the evolution timescales for stars in different stages. Finally in ^we discuss 
the body of evidence involving second parameter effects (cluster to cluster variations) and 
intracluster differences between stars. 



2. Datasets 

2.1. Archival Hubble Space Telescope Imaging 

We used WFPC2 images taken in three different studies and ACS images from one 
additional study for photometric measurements. The principal investigators and filters used 
are listed in Table |2j The images from proposal 8174 were composed of three partially 
overlapping fields that were reduced separately and later combined. 

All of the WFPC2 i mages wer e proces sed using the HSTPhot@ package (version 1.1.7b), 



which was described by iDolphinI fcOOOaf l. Star positions for each set of processed data 
were derived using the IRAF task METRIC, which we used to convert the HSTPhot pixel 
coordinates to RA and DEC. METRIC positions often have absolute errors of up to O'.'S, 
but absolute errors are not important for our purposes. The RA and DEC coordinates were 
converted to a system relative to the cluster center (Fig. [T]). The data for each star, including 
coordinates, flight system magnitudes, and errors, were extracted from each dataset and 
merged into one master file. In cases where a star had multiple measurements in the same 
filter, the values were averaged using weights derived from the measurement errors. This 
master dataset was used to construct CMDs that we used to identify post-MS stars. The 
F160BW images were reduced separately because the filter distorts star positions relative to 
the rest of the images, complicating the matching of stars within HSTPhot. After completing 
the photometry, we were then able to match all of the stars with previously identified HB 
and AGB stars. 



^http : //purcell . as . arizona. edu/hstphot/| 
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The ACS Wide Field Camera (WFC) images were processed using the DOLPHOli 
package with its module for ACS data. DOLPHOT tasks mask bad pixels, correct for effective 
pixel area, and do point-spread function photometry. Individual frames (prior to drizzling in 
the ACS pipeline) were obtained from the HST Archive in order to get photometry on stars 
covering the widest possible dynamical range (including stars that were saturated on long 
exposures) and to avoid the pixel resampling that goes along with the drizzling process. 

Because the WFC sits far from the optical axis of HST, there is significant geometric 
distortion of these images. For the purposes of the photometry, the differences in effective 
pixel area were corrected within DOLPHOT through the use of a pixel area map provided 
on the ACS websit^. DOLPHOT, however, does not currently correct pixel positions to sky 
coordinates. For the purposes of relative astrometry, we used the fourth-order polynomial 
corrections provided in the most recent IDCTAB file for the dataset (qbul641sj Jdc.fits). 
These corrections reduce distortions to about 0.1 pixels, which is more than adequate for our 
purposes. The star positions were later put on a common coordinate system with photometry 
from WFPC2 images as described below. 



2.2. Archival CFHT Images 

The Hubble Space Telescope is very effective at resolving the core of the cluster, but 
its approximately 3' by 3' field of view is insufficient to capture most of the cluster stars. 
We thus turned to the archive of the Canada-France-Hawaii Telescope (CFHT) CFH12K 
camera for data covering the majority of the globular cluster (Fig. [1]), with data collected 
by J.-H. Park, Y.-J. Sohn, and S. J. Oh in Feb. 2001 using B, V, and Cousins / filters. This 
camera is composed of a mosaic of 12 CCD's, of which we only used numbers 11 and 12 
because these two chips covered the most heavily populated parts of the cluster. Each chip 
covers a roughly 7' by 14' area of sky. The images we obtained from the CFHT archive 
were already bias-subtra cted and fiat- fielded. For the photometry, we used the DAOPHOT 



II/ALLSTAR packages (IStetsonI 119871 ). Star positions in the CFHT images, like the HST 
images, needed to be transformed from pixel coordinates into angular offsets from the cluster 
center. We cross-identified stars from the USNO A2.0 catalog, covering an approximately 
18.3' by 16.7' area of sky around the cluster that encompassed the CFHT field. These 
catalog stars were used to convert the CFHT pixel positions into RA and DEC coordinates 
relative to the core of M13. In regions of overlap, the HST datasets were placed on the same 



•^ ihttp : //purcell . as . arizona.edu/ dolphot/ 
http : //www. stsci . edu/hst/acs/analysis/PAMS 
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coordinate system. 



2.3. KPNO 0.9m Imaging 

We also reduced imaging the KPNO 0.9 m telescope taken in B, V, and Cousins / 
filters by Bolte and Sandquist on two nights in May 1995. The camera used a single 2048 x 
2048 CCD with a scale of 0'.'68 per pixel, for a field of view 23'2 on a side. The photometry 
was done similarly to the CFHT dataset, and coordinates were transformed to the CFHT 
system. 



2.4. Proper Motion (PM) Studies 



We cross identifie d all stars from the proper motion studies of lCudworth &: Monetl (119791 ) 
and ICudworthI (119791 ) . The field covered by these studies is similar in size to that of the 
KPNO images. For the brightest stars, membership probabilities were available to within 
about 1' of the cluster center, while for fainter stars the innermost 2 — 3' was effectively 
excluded. Their photometry had a faint limit at V ^ 15.6, which is about 0.5 mag fainter 
than the red end of the HB. As a result, the proper motion membership probabilities were 
useful primarily for identifying bright stars that would otherwise contaminate post-HB and 
giant star samples. 



2.5. Archival UIT Imaging 



The archive of the Ultraviolet Imaging Telescope contains aperture and point-spread 
functi on photometry of M13 produced following the procedures described in IStecher et al. 
(Il997l ). and the images involved have previously been discussed by lParise et al.l (119941 ). How- 
ev er, in examini ng the UIT images, we found a large number of sources that were not recorded 
by iParise et al.l or in the UIT archive file, and these match up with horizontal branch stars 
identified with optical photometry. It appears that the earlier photometry had an unneces- 
sarily high faint limit. 

We therefore rereduced the two longest exposure far UV images (using the B5 filter with 
central wavelength 1620 A) taken on the Astro 2 mission (fuv2418, 192.5 s; fuv2419, 953.5 
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s). We downloaded archived IDL reduction routines called MOUSSElf] for use on linearized, 
fiat-fielded digitizations of the original photographic images. The uit_find routine identifies 
significant sources using an algorithm based on DAOPHOT, and determines centroided posi- 
tions for each. The uit_ape r routine coiiducts aperture photometry. We used a 3 pixel radius 
aperture recommended by iParise et al.l (119941 ) and verified that the measurements matched 
the archived aperture photometry to within 0.01 mag for sources in common. The shorter 
fuv2418 exposure was analyzed to allow the brightest sources to be incorporated, as they 
were saturated on the longer exposure. 

Many sources could be identified in the core of the cluster because crowding is relatively 
mild in the far UV. However, crowding does significantly modify the photometry, so we 
will generally re s trict o ur discussion to stars more than 120" from the cluster center, as 



did iParise et al.l ( 119941 ). This partly resulted from the pixel scale (pixels in the digitized 
photographs have a scale of about 1'.'14) and tracking inaccuracy in the longer fuv2419 
exposure. We identified a few instances of blending in the UIT observations outside this 
radius when stars on the blue tail of the HB were close enough together on the sky. We 
rated UIT detections by hand based on the apparent degree of contamination, and will only 
plot stars with the best ratings. 

We note that only in one case (ID 72) did we find a star that was very likely to be a 
proper motion nonmember, and in two other cases (UIT source 568; ID 3904) did we find 
stars that could not be identified as HB or AGB stars. ID 72 has optical photometry that 
places it only about 0.2 mag brighter than the HB level. We found no clear counterpart to 
UIT source 568 (mi62o = 15.98), an d hypothesize t hat it might be related to the faint UV 
sources found in HST photometry by lFerraro et al.l (119971 ). Source (ID 3904) corresponds to 
a bright blue straggler star in optical photometry. 



2.6. Photometry Calibration 

2.6.1. CFHT 



We calibrated the BVI PSF ph otometry from CFHT to the standard system using 
standard stars from M13 described in IStetsonI (120001 ). These standard stars are rigorously 
matched to the same photometric system as the earlier iLandoltl (119921 ) study. We obtained 
the 2 May 2005 update star list from the CADC websit^. The conditions for the CFHT 



■"^http: / /archive. stsci.edu.uit / analysis.html 
''http : //cadcwww . dao . nrc . ca/stEindards7 
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observations did not appear to be photometric, so the magnitude zeropoints for each image 
were free parameters in the cahbration. 

Standard stars in the Stetson field cover almost the entire range of colors for evolved 
stars: from stars at the faint end of the HB distribution to a giant star within 0.5 V mag of the 
RGB tip. Three redder stars (probably field) were included to extend the color calibration 
closer to the color of the tip {B — I ^ 3.1). We found it necessary to calibrate the two 
CFHT12k CCD samples separately because of significant differen ces in color terms. (This is 



consistent with the results of oth er studies — for an example, see lKalirai et al.ll200ll .) Using 



the CCDSTD (e.g.. IStetsonlll992l ) program, the standard star transformation equations were 
found to be: 

b = B + ai + (-0.025 ± 0.006)(5 - /) + (0.0084 ± 0.0029)(S - if 

v = V + bi + (-0.0038 d= 0.0026)(fi - /) 
i = I + Ci + (-0.0079 ± 0.0025)(fi - /) 

for chip 11, and 

b = B + ai + (-0.030 ± 0.004)(5 - /) + (-0.0173 ± 0.0018)(5 - if 

v = V + bi + (-0.0084 ± 0.0013)(5 - /) 

i = I + Ci + (-0.0011 ± 0.0012)(5 - /) 

for chip 12, where b, v, and i are instrumental magnitudes, B, V, and / are standard 
magnitudes, and Oj, 6j, and Cj are the zero points for individual frames. The {B — I) color 
was chosen primarily for its wide wavelength spacing, which helps minimize the importance 
of photometric errors in B or /. We compared the final calibrated measurements with the 
Stetson standard values, as shown in Fig. |2J The comparisons include 417 stars from chip 
11 and 286 stars from chip 12. The median residuals AB, AV, AI, and A{B — I) were all 
less than 0.003 mag for the samples from both CCDs, and consistent with zero. There do 
not seem to be significant systematic trends at the extremes of the color range either. 



2.6.2. HST 

For most purposes, we preferred to leave the HST datasets in flight system magnitudes 
using the VEGAMAG zeropoints in order to maintain the relative precision of the original 
photometry and to make it possible for others to use improved transformations to a stan- 
dard system in the future. However, to produce uniform samples of stars from photometric 
datasets in different regions of the cluster, we did need to derive some transformations. 
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Because the CFHT photometry appears to be accurately matched to the standard sys- 
tem and because the CFHT images had good spatial resolution in the cluster core, we use this 
dataset to examine the calibration of the ACS WFC photometry in V and /. The comparison 
involves HB and faint RGB stars because stars high on the RGB were heavily saturated even 
on the shortest exposures. In comparing the ACS flight system magnitudes to the CFHT 
photometry, it was clear that second-order color terms in color were needed to transform 
the F606W observations, but a first order term seemed sufficient to transform the F814W 



observations. This is in agreement with the t ransformations derived by lSirianni et al.l ( l2005l ). 



However, we found that when we used the Sirianni et al. transforms, the RGB stars were 



systematically too faint by about 0.04 mag. We therefore derived our own transformations. 
It was necessary to iterate the process, eliminating stars with large positive residuals (most 
likely due to blending in the CFHT dataset). The following transformations were derived 
from 437 stars: 

meoe = V + (0.010 ± 0.004) + (-0.0592 ± 0.0174)(F - /) + (-0.1553 ± 0.0207)(r - I)^ 

msu = 1 + (-0.005 ± 0.004) + (0.0290 ± 0.0053)(V - /) 

Fig. [3] shows the comparison between the CFHT photometry and the ACS photometry 
calibrated to VI. 

For a small number of stars, WFPC2 images were the primary source for photometry, 
thanks to crowding and/or gaps in sky coverage. We used WFPC2 photometry in the 
VEGAMAG system, but have calibrated measurements in the F555W an d F785LP bands 



to V and /. We first used the standard calibration of iDolphinI (l2000bl ). but found that 
there were large color-dependent residuals. As a result, we determined new transformation 
equations using 153 RGB and HB stars in common with the CFHT data (eliminating stars 
that were most likely to be blended, as we did with the ACS data): 

^555 = V + (0.065 ± 0.009) + (0.0860 ± 0.0364)(\/ - /) + (0.0207 ± 0.0498)(l^ - if 

^785 = / + (0.051 ± 0.010) + (0.1346 ± 0.0413)(y - /) + (-0.2869 ± 0.0557)(V - J)^ 

The transformed WFPC2 and CFHT datasets are compared in Fig. |H Although the formal 
errors on the quadratic term for the msss equation imply it has low significance, it was our 
qualitative judgement that it improved the transformation at the red end of the color range. 

For a small number of stars at the faint end of the RGB in the cluster core, the only 
reliable near-infrared measurements were in the WFPC2 field of proposal 8278, and so this 
was also calibrated to CFHT measurements using a sample of 87 RGB and HB stars: 



^555 = V+ (0.052 ± 0.011) + (0.0538 ± 0.0428)(y - /) + (-0.0378 ± 0.0629)(\/ - if 
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msM = / + (0.021 ± 0.012) + (0.3233 ± 0.0481)(l^ - I) + (-0.3328 ± 0M8S){V - if 

In the color range covered by the stars with WFPC2 photometry, the two cahbrations typi- 
cally agree to better than 0.05 mag. 



2.6.3. KPNO 

We observed Landolt standard fields on one photometric night from the KPNO run, 
and a standard s tar calibration using these data was used in calibrating cl uster data from 



previous studies (IHargis et al.ll2004i iPoUard et al.ll2005l : iFekadu et al.l 120071 ). In this paper 



we have opted to calibrate the KPNO photometry directly against the well-observed Stetson 
(2000) standards in the M13 field used to calibrate the CFHT data. To ensure the calibration 
of the evolved populations, we restricted the standard star sample to 343 stars brighter than 
the approximate turnoff of the cluster {V < 19). The transformation equations are 

b = B + ai + (-0.0279 ± 0.0018)(5 - I) 

v = V + bi + (-0.0002 ± 0.0017)(S - /) 
i = J + c, + (-0.0130 ± 0.0020)(5 - /) 

The residuals of the calibration are shown in Fig. [51 The median residuals in all cases 
{B, V, I, B — V, and V — I) are less than 0.005 mag, and well within the uncertainty in the 
fit. There may be some small systematic trends at the extreme ends of the color range for the 
standard stars: {B — V) < —0.1 and (-B — V^) > 1.2. After additional experiments though, 
we found that second-order color terms or changes to the standard star sample (fewer main 
sequence turnoff standard stars, so that the giant and horizontal branch stars carried more 
weight) did not significantly improve the fits. 

For most of the scientific applications later in the paper, the relative precision of photo- 
metric measurements within one dataset is most important. However, for the discussion of 
HB and RGB luminosity functions, it is important that there are not large systematic errors 
between the datasets we have merged together. The luminosity functions will employ the 
J-band measurements, so particular attention should be focused on that filter. Figures [6] and 
[7] show comparisons between the KPNO photometry and the CFHT photometry from chip 
11 and chip 12. There appear to be signs of systematic residuals in various places [notably 
some small offsets in V and / in chip 11, and a trend in A (5 — V) versus {B — V) for chip 
12]. However, for the great majority of stars, the residuals are quite small (< 0.04 mag). In 
the KPNO-CFHT comparisons, the /-band measurements seem to be well calibrated. 
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3. Identification of Evolutionary Status 

Our first step was to fi r st ide ntify all known variable stars (RR Lyrae, BL Her, bright 



giants) from iKopacki et al.l ( 120031 ) in order to eliminate stars that might confuse classifica- 
tions. Our subsequent procedure for separating stars by evolutionary stage involved all of 
the datasets described above (from widest area coverage to smallest: PM, KPNO, CFHT, 
HST ACS, and HST WFPC2). We relied primarily on the observations with the highest 
available spatial resolution for each portion of the field. 

The main CMDs used to identify giant stars (RGB and AGB) are shown in Fig. [HI For 
the outermost portion of the field, the KPNO dataset was used along with proper motion 
information. The CFHT dataset was used for the bright RGB and AGB samples even within 
the core, driven by the fact that the brightest cluster stars were saturated in the HST images 
(with the exception of some exposures in the ultraviolet). Crowding effects were minor for 
the CFHT images, which were taken with good seeing and with a high resolution camera 
(0'.'2 pixels). Because giant stars (RGB and AGB) were all observed in BVI filters and 
because the B — I color covers the widest wavelength range, we used B — I as their primary 
discriminant. For fainter RGB stars (J > 14.5), HST photometry was available in the core, 
and so ACS WFC photometry was used when available, and WFPC2 photometry for all 
other stars falling between WFC chips or outside the WFC field. 

The HB and AGB manque are comprised of hotter stars, so we relied on the shortest 
wavelength filters to identify them. For the portions of the core observed in the ultraviolet, 
the F160BW, F255W, and F336W filters on HST are particularly good discriminants (see 
Fig. [H]), with the F160BW observations effectively selecting only the hottest HB stars and a 
handful of manque stars. Otherwise the HST ACS observations were the primary source for 
identifications for the bluest horizontal branch stars in the core. HB candidates in the core 
were selected based on CMD position, but derived x^, sharpness, and crowding values were 
used to fiag candidates that might be spurious. Questionable candidates were also examined 
by eye, and stars that clearly fell on chip defects or very close to brighter stars were rejected. 
The crowding parameter in DOLPHOT (the change in measured brightness if close neighbors 
had not been simultaneously fit with point-spread functions) was very effective at flagging 
spurious stars. We found that a change of 0.4 mag indicated a spurious star in nearly all 
cases. The optical CMDs used in the HB star selections are shown in Fig. (TDl 

Just outside the observed HST flelds, there is noticeable scatter in the photometry 
measured from the CFHT images, as can be seen in the lower right panel of Fig. [TOl 
To reliably identify blue HB stars in this region, we made use of UIT observations. UIT 
observations cover almost the entire fleld discussed in this paper, and even though blending 
sometimes made accurate photometry a problem in the core, the UIT observations allowed 
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us to identify blue HB stars even there. As a result, we believe we have virtually all of the 
faint HB stars in our sample, all the way from the cluster center out to the edge of our large 
ground-based fields. 

AGB manque stars also deserve discussion. In optical CMDs these stars (post-HB stars 
that will not reach the traditional AGB) are predicted to evolve almost parallel to the zero- 
age HB (ZAHB), and the hottest ones can sometimes be hard to distinguish from normal 
HB stars. Ultraviolet-optical CMDs can flatten out the blue HB and make it possible to 
identify AGB manque stars by their relative brightness. For the majority of the field we 
identified candidates from the HB sample in the highest resolution optical CMD first. We 
subsequently used observations in F160BW and F255W for the core, and the wide-field 
UIT observations elsewhere to validate most of the candidates. We found that nearly all 
of the bluest candidates selected from the optical CMDs were verified with the ultraviolet 
observations, but 5 additional manque stars were identified first in the UV. For the UV 
observations, we defined stars to be manque if they were more than about 0.8 mag brighter 
than the faintest HB stars at the same color. (This criterion is based on where theoretical 
tracks predict evolution slows again after core He exhaustion.) The UIT observations were 
subject to blending in the central regions and some candidates were rejected based on high- 
resolution optical photometry, but in all cases it was possible to identify the stars that were 



the pr imary UV sources based on other photometry. In addition, some of the IWhitney et al. 



(Il995l ) sources did not appear to have cluster counterparts at the tabulated positions. We do 
not have an explanation for this. Two stars (HB 72 and 295) present reasonably strong cases 
in optical CMDs for being manque, but are not discernibly unusual in the UV, so we have 
left them as HB stars. Finally, two manque stars (AGB 12 and 26) were identified using 
only optical photometry because they fell outside the HST ultraviolet observations, but in a 
region of the cluster where the UIT images were too crowded to yield reliable measurements. 
All of our cross identifications are given in Table O The photometry for the identified stars 
can be found in Tables HHIl With the use of the UV observations, we believe we have isolated 
nearly all AGB manque stars from the HB sample. 



4. The Red Giant Branch 

We assembled a comprehensive tabulation of the bright RGB stars in M13 from the tip 
of the RGB to the base of the RGB, where it meets the subgiant branch. As described in 
^ we eliminated AGB stars from the sample using the best photometry available (see Fig. 
[8]). This is a relatively easy task in M13 because the AGB is well-separated in color from 
the RGB except near the bright end. 
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In a study of the cluster core using the HST Planetary Camera, ICohen et al.l (119971 ) 
previously found evidence that there was an excess of bright giants compared to models. 
Though their photometry was not precise enough to separate AGB and RGB stars, the excess 
appeared to be too large to be explained by the AGB star population. To examine the evo- 
l ution rates of the bright giants, we followed the procedure described by ISandquist fc Martel 
( 120071 ). In this method theoretical cumulative luminosity functions in / band are shifted 
in magnitude so that the tip of the giant branch matches the brightest giant, and the lu- 
minosity fun ction is then normal ized to the total number of stars found brighter than the 
RGB bump. ISandquist fc Martell found that this has the benefit of making the comparison 
extremely insensitive to chemical composition and age, which makes it possible to test the 
physics of the upper giant branch in clusters with large giant populations. In that study, 
Sandquist &: Martell found that the giant population in NGC 2808 was significantly depleted 
relative to the model predictions, and hypothesized that this could be due to underestimated 
neutrino emission rates or loss of the giant envelope (causing some giants to leave the RGB 
before undergoing a normal helium flash). We have therefore tested w hether M13's very 
blue H B morphology is reflected in the character of the giant population. ISandquist fc Hess 
( I2OO8I ) recently examined a large sample in the core of the distant cluster NGC 2419 (another 
cluster with a long blue HB tail), which showed only weak signs of a depleted bright giant 
population. 



To c onduct the test for M13, we shifted Victoria- Regina (IVandenBerg et al.ll2006l ) and 
Teramo (ICassisi et al.l 12004 ) models to the / magnitude of the brightest giant. We need to 
be careful because almost all of t he brightest giants in M13 are known to vary irregularly 
on timescales of more than 40 d (IKopacki et al.l 120031 ). The variabi lity itself shou l d be a 
relatively small effect based on the variation amplitudes observed by iKopacki et al.l ( l2003l ): 
0.07 < AV < 0.38 for the 12 stars studied, with most having amplitudes toward the low end 
of that range. We identified all of the known RGB variables, and find that VI 1 appears to be 
the brightest (/ ~ 10.3) in our study and most others. (Our l^-band measu rements appear 



to have been made in brighter periods in the star's variation. According to IKopacki et al. 



the star's mean V magnitude is only a few hundredths of a magnitude different from other 
bright giants.) 

Fig. im shows the comparison for a model having [Fe/H] = —1.41 and [a/Fe] = -1-0.3. 
We tested models with different metallicity, but as described earlier, the exact value does not 
significantly affect the results. In addition, the Victoria-Regina and Teramo models agree 
well with each other. The plotted Victoria-Regina model shows the best agreement with the 
observed position of the RGB bump. We conducted Kolmogorov-Smirnov (K-S) tests with 
the observed and theoretical luminosity functions using different faint cutoffs for the RGB 
sample. (Because the RGB sample increases exponentially with increasing magnitude, this 
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was necessary to ensure that deviations on the less-populated bright RGB are not washed 
out.) The test results are presented in Table O The probability P that the observed stars are 
consistent with being drawn from the theoretical distribution reaches a minimum of about 
5% when the cutoff is at / = 12.76 in response to a clump of stars at / ~ 12 (see below). 
Otherwise, the probabilities remain greater than 15%, and the overall impression from the 
test is that the re is at most aslight deficit of stars compared to the model predictions, in 
contrast to the ICohen et al.l (119971 ) results. 



As discussed by ISandquist fc Martell ( 120071 ) . this method of comparing with models 
is insensitive to chemical composition for [Fe/H] < —1, which means that it should be 
possible to combine samples from different clusters to improve the sta tistical significance 
of th e comparisons. M5 has a readily available sample of RGB stars ( ISandquist &: Bolte 



20041 ). and issues related to combining the two samples are minimal. For exa mple, M5 and 



M13 have ver y similar distance moduli [(m — M) 



1999; -0.02, 



Recio-Blanco et al. 



A/5 



-0.01, 



Harrii 



(m 



M) 



0.06, Ferraro et al. 



2005 



0.08, ICarretta et al.ll2000i and r eddening [E(B 



1996; -0.02, 



Recio-Blanco et al. 2005 values, and their 



V)Mn-HB-V) 

metallicities are similar (with M13 universally identified as more metal poor by about 0.2 
dex). The M5 sample is similar in size to the M13 sample as well. In the comparisons below 
we use a model with a compromise [Fe/H] = —1.31, although once again the composition 
choice is not critical. In light of the similar and relatively well-determined distance moduli, 
we have chosen to combine the samples after correcting for the small difference in distance 
modulus. This is probably superior to shifting the cluster samples so that the brightest 
giants in each cluster match because i) many of the brightest giants in M13 are variable, and 
ii) statistically speaking, the brightest giants may fall at different brightness levels fainter 
than the tip of the red giant branch (TRGB). 

The combined sample is shown in Fig. [T21 K-S test probabilities only go below 10% 
in small ranges of magnitude (for example, P = 0.05 for a cutoff at / — Itbgb = 3.00). 
Overall, we regard the agreement with t he models as good, and a confirni ation of the most 
recent plasmon neutrino emission rates ( iHaft et al.l Il994j : lltoh et al.l Il996l ) , which influence 
the evolution rates near the TRGB. There are variations in counts around the theoretical 
predictions, but are not significant at more than a 2a level. With these samples, we see no 
reason to consider particles beyond the standard model with very large mean fr ee paths (such 



as ax ions and WIMPS) that could affect cooling of the degenerate core (e.g. ICatelan et al. 



19961). 



The evolution rates for fainter giants can al so be exarn i ned u sing star numbers. Using a 
sample covering a smaller portion of the cluster, ICho et al.l ( 120051 ) found that M13 appeared 
to have "extra" stars in and slightly fainter than the RGB bump. On the other hand, in 
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their study of the bumps of a large sample of clusters, iRiello et al.l ( 120031 ) found that M13 
appeared to have a lower than expected number of bump stars (stars with V within ±0.4 
mag of the bump center ) compared to gi ant stars between .5 and 1.5 V magnitudes fainter 
than the bump. Both iRiello et al.l and iBono et al.l ( 120011 ) found Ml 3 had a lo wer value 
than the mean for their whole sample of clusters, and in the case of IRiello et al.l the value 
was almost 3cr lower. Usi ng Vbnm.p = 14.75, our M13 samples are more than twice as large 
as the ones tabulated bv IRiello etaP We &nd Rh„mp = 29 7/609 = 0.487 ± 0.034, which is 



uite consistent wi t h the mean value found by IRiello et al.l and with theoretical predictions 
IRiello et al.l l2003l : iBjork fc Chaboyerl |2006[ ) . The slope of the lower giant branch in the 
cumulative luminosity function (which is related to the evolution timescale) for M13 is also 
in very good agreement with the predictions of the Victoria-Regina and Teramo models. 
Based on our large sample, the M13 giants are evolving at the rate predicted by models. 

As mentioned above, the CLF of M13 also shows a slight enhancement in the number 
of RGB stars at J ~ 12.0 ot V ^ 13.15. This feature can be seen by eye in some of the 
CMDs ( for example, see F ig. [T3|) . This is reminiscent of the "heap" feature seen in NGC 
2808 by lBedin et al.l (120001 ). which was identified about 1.4 magnitudes brighter than the 
RGB bump. This feature falls approximately where the RGB and AGB begin to overlap. 
While there have been quite a number of spectroscopic studies of M13, there is not an 
unambiguous spectros copic signature to separate AGB stars from RGB stars. For example. 
Smith &: Brileyl ( 120061 ) find that photometrically-identified AGB stars had uniformly weak 
CN bands, indicative of O — >N conversion. However, there is a significant group of CN- 
weak RGB stars even fainter than the level where they could reasonably be confused for 
AGB stars. 

We present the differential luminosity function for the bright end of the RGB in Fig. [T4l 
As iBedin et al.l found for NGC 2808, the heap is not statistically significant. Its appearance 
in the CMD comes largely from the contrast with slightly fainter RGB stars, which appear 
to be less common than predicted theoretically. In our M13 sample (which is close to as 
large as it will get for this cluster), the heap appears to be significant at about the 2.5o" level 
based on K-S tests on the cumulative luminosity function. 

To summarize, after gathering the largest possible sample of giant stars in the cluster, we 
see at most low-significance deviations from theoretical luminosity functions. We emphasize 
that our luminosity function comparisons have been done in a relative sense, by forcing 
agreement between the brightness of the observed and theoretical TRGB. In contrast to the 
case for NGC 2808 (but in agreement with NGC 2419), we do not see evidence of large 
numbers of stars leaving the RGB early that could account for the bluest HB stars in the 
cluster. 
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The Horizontal Branch 



5.1. Brightness and Radial Distributions 

M13 possesses a horizontal branch comprised mainly of stars bluer than the instability 
strip (Fig. [TOj) . Before we turn to comparisons of the HB population with other post-main- 
sequence samples, we will first examine the distribution of HB stars in the CMD and as a 
function of radius in the cluster. 

Our dataset has been derived from a wide-range of source material, but V photometry 
is common to all of the so urces, and / photometry is available for all but one of the sources 
(the proper motion data of lCudworth fc Monetlll979l ). We use I photometry to describe the 
relative positions of stars on the HB because it is an observable quantity, and HB stars grow 
monotonically fainter in I with decreasing mass. The relationship between / magnitude and 
mass results from competition between decreasing radius and increasing surface temperature 
as the mass of the star's envelope decreases, and from the rapidly increasing bolometric 
corrections that result. The use of a directly observable coo rdinate for HB position also 



reduc es the difficulties in comparing models and observations (IVargas Alvarez fc Sandquist 
2007h . 



Because of the potential influence of blending, we drew the / magnitude of a star from 
the photometric dataset with the highest spatial resolution: ACS WFC, CFHT, and KPNO. 
Five stars fell in the gap between the two ACS WFC chips, and their photometry was 
derived from WFPC2 observations in F785LP ( §2.61) . In Fig. [TBI we show a histogram of 



the HB star distribution in /, clearly showing the well-known bimodality. For conceptual 
purposes, we have broken the HB into three parts: / < 16.25 (the bright peak, hereafter 
PI), 16.25 < / < 18 (intermediate stars), and / > 18 (the faint peak, hereafter P2). The 
break at / = 18 corresponds to gap G3 identified by lFerraro et al.l (Il998l ). while the break at 
/ = 16.25 is slightly fai nter than their g a p Gl a nd at approxima t ely th e same position as the 
"m jump" identified by iGrundahl et al.l (119981 ). iGrundahl et al.l ( 1l999l ) tentatively identified 
this as the place on the HB where radiative levitation of heavy elements becorn es important, 
and spectroscopic studies (IBehii l2003l : iFabbian et al.l l2005l : iPace et al.l l2006l ) confirmed a 
change in the atmospheric abundances of heavy elements at similar positions in a number of 
other clusters. The gap Gl appears to correspond to a loca l minimum in the HB distribution 
at / ^ 15.7. There is no doubt that iFerraro et al.l (119981 ) saw clear evidence of gap Gl in 
their dataset, but the gap in their data was fairly wide and it is difficult to assign a precise 
location using their smaller dataset. 



Because the well- studied cluster MS s hows a clear change in HB morphology with dis- 
tance from the center (ICatelan et al.ll200ll ). we examined the distribution for signs of radial 
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variations. In the top panel of Fig. [151 we see that the HB and RGB distributions are essen- 
tially identical. In the lower panel, we plot the HB distributions for samples with r < and 
r > (380 and 404 stars, respectively). There appear to be slight differences between the 
sample distributions, with the second peak at the end of the blue tail becoming less promi- 
nent in the outer portion of the cluster. As shown in the cumulative radial distributions 
in the bottom panel of Fig. [161 there are modest differences in the innermost 140" of the 
cluster, with the intermediate HB stars being the least centrally concentrated and the faint 
peak stars being the most centrally concentrated. However, for this global comparison, a 
Kolmogorov-Smirnov test indicates that there is still a 21% chance that those two distribu- 
tions could be drawn from the same parent population. When the comparison is restricted 
to r < Tc or r/i, the probability does not change significantly. 

Table [H] shows the number counts for different radial samples. The counts show a 
tendency for stars in the faint peak to be more centrally concentrated. They are almost as 
abundant as bright peak stars in the very center of the cluster (r < rc/2), but the fraction 
drops quickly so that for r > Vc, they are consistently about as abundant as the intermediate 
group. There is a marginal trend in the fraction of HB stars in the faint blue peak (/p2) as a 
function of distance, but a significant change in the difference in the fractions in the bright 
red and faint blue peaks (/pi — fp2)- 



Cohen et al.l (119971 ) found that the blue HB stars in M13 appear to be centrally depleted 
relative to other types of stars for their sample with r < 60". We believe the difference is 
probably due to incompleteness at the faint end of their horizontal branch sample. The 
appearance of the horizontal branch in their Fig. 4 clearly shows the bright peak, but there 
is little or no sign of the faint peak (most stars in 18 < V < 19), which falls near the 
completeness limit of their dataset. 



5.2. Notable Stars 

Only nine RR Lyrae stars are known in the entire cluster, but there are stars redder 
than the instability strip. We find a conspicuous group of eight stars just bluer than the 
RGB but approximately 0.5 mag fainter than the AGB clump (seven are seen in Fig. [TT]) . 
A star at the red end of the theoretical ZAHB locus plotted in Fig. [TTjhas a mass of O.8M0, 
so HB stars in this part of the CMD would be more massive than stars at the cluster turnoff 
(the HB is theoretically expected to curve back toward higher temperature and luminosity 
for more massive stars). We believe that th ese stars are probably e volved counterparts to 



blue straggler stars, as originally proposed by lFusi Pecci et al.l (119921 ). This possibility could 



be tested by comparing the number to the size of the sample of cluster stragglers (and 
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thereby comparing their relative lifetimes), but this is beyond the scope of the current work. 
Although field stars are relatively common in this part of the CMD, five of the eight star s 



have proper motions indicating that they are high probability members (ICudworthl Il979l ) . 
The remaining three stars are all within 45" of the cluster center. Only one of the stars (HB 
532) has UV observations — in the F255W passband it is significantly brighter than the 
RGB but fainter than hotter HB stars, which rules out an optical blend of a hot HB and an 
RGB star. 

From careful examination of the photometry, we identified three additional stars that 
we believe may be red HB stars. HB 314 has proper motion information that identifies it as 
a high probability member and photometry that places it at V magnitude nearly identical to 
stars blueward of the instability strip. The two remaining stars do not have proper motion 
information, and their CMD positions give them greater probabilities of being field stars. HB 
793 has a V magnitude placing it at the cluster HB level, and it is fairly blue compared to the 
known field stars of similar brightness, but it is outside the proper motion field (r = 625"). 
HB 133 is closer to the cluster center (r = 103"), but has a V magnitude placing it about 
0.3 mag fainter than the cluster's HB. However, there are no other stars near it in the CMD, 
and there is a small possibility tha t it is an undetected variable star, although it was within 



the field studied by iKopacki et al.l (|2003[ ) using a sensitive image subtraction method. 



Although quite a few globular clusters have long blue horizontal branch tails (NGC 
6752, M2, and M80 are examples), M13 is somewhat unusual in having a secondary peak 
in its distribution of stars that falls near the theoretically-predicted blue end of the HB. 
Clusters that have blue hook stars (such as u Cen, NGC 2808, and NGC 2419) have blue 
HB tails in optical CMDs that extend to fainter magnitudes. As a result, we have looked at 
the faintest HB stars in the optical to gauge whether there are any blue hook stars present. 



Sandquist fc Hesd (120081 ) noted that the M13 HB largely terminates at the approximate 
position of a gap in the EHB population of NGC 2419 when the HBs are aligned at their 
brightest points in B. In CMDs employing B, F439W, or F336W filters, we identified a 
handful of candidates that were fainter than the faint HB peak P2. However, ultraviolet 
CMDs must be used in conjunction with the optical CMDs in order to get a full picture 
of the spectral output of the stars. In Fig. [HI we present our UV CMDs for M13. As 
we progress from the near UV (F336W and F255W) to the far UV (B5 and F160BW), 
bolometric corrections become less severe for the hottest stars, and the hottest HB stars 
match the rest of the blue HB in brightness. Although the WFPC2 F160BW and UIT B5 
filters cover similar portions of the spectrum, and produce similar CMDs in Fig. [HI the 
F160BW filter reaches peak sensitivity more than 100 A shorter than the B5 filter. We 
believe this is responsible for the slight upturn near the blue end of the HB in F160BW. In 
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CMDs using even shorter wavelength observations, the hottest HB stars would continue to 
change their configuration. 



D'Cruz et al.l (|2000[ ) found a large group of blue hook stars in uj Cen that were sublu- 
minous in an (F160BW, F160BW-V^) CMD compared to other stars at the end of the HB. 
A similar structure does not exist in Fig. [9] near F160BW— F555W~ —3.5, but there are 
two stars that meet the typical definition of the blue hook: fainter than the ZAHB in the far 
UV. The optical measurements for these two stars put them significantly fainter than the 
predicted hot end of the HB. The UIT source (HB 2) is approximately 0.7 mag fainter than 
the hot end of the HB in the B5 filter, and 1 magnitude fainter in B. (It is the faintest HB 
source in the KPNO panel of Fig. [TOl) HB 431 is among the faintest sources in the three 
WFPC2 UV filters, but is faintest only in F160BW. (It is the second faintest HB source in 
the ACS WFC panel of Fig. m) 

Nine other stars are fainter than the EHB in optical bands and sit at the extreme 
blue/faint end of the HB in ultraviolet bands (although fairly close to the extrapolated zero- 
age HB), and so we fiag them as potentially interesting stars. The UIT sources are HB 8, 
56, 73, 86, and 757, while the WFPC2 sources are HB 373, 485, 611, and 636. 



5.3. Evolution of HB Stars 



With a large sample of HB stars, the sample of evolved HB stars also grows, making it 
possible to trace post-HB evolution. The two most densely populated portions of the HB 
contain 75% of the stars, and so to first order, post-HB stars will trace the evolution of these 
two groups. 



In Fig. [T8]and[T9l we show a comparison betwe en theoretical models of lPietrinferni et al. 



(120061, hereafter, Teramo) and iDotter et al.l ( 120071 . hereafter, DSEP) and our ground-based 
photometry in (V, i? — /) for r > 200". To gauge the edges of the HB distribution, we used 
zero-age and central helium depletion {Y(. = 0.05 or 0.10) CMD positions as a function of 
mass. Because there is som e uncertainty about the precise chemical composition of M13 



(e.g., iJohnson fc Boltd Il998l ). we fitted the HB band to the photometric data in order to 
better identify evolved stars in a differential sense. We emphasize, however, that the fit does 
not completely validate the models. In fact, there are mismatches between observation and 
theory that we will examine below. 
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5.3.1. The Reddest BHB Stars 



Because M13 has a blue HB, traditional photometric indicators of distance {My ), age 
{AVr^(f, for example) and helium {R = NhbINrgb) that have the V magnitude of the HB 
in their definitions cannot be used without corrections of uncertain accuracy. As a result, 
redder stars may better define the brightness of the HB if they can be proven to reside near 
the ZAHB in the CMD. 



As shown in Fig. |20l we identified 55 stars redward of the most heavily populated part 
of the HB (the bright peak PI; 308 stars in 16.25 > / ^ 14.9). The stars stretch to the 
instability strip (0.10 < V — I < 0.22), and it is likely that at least some of the RR Lyrae 
stars belong to the same group. These stars have largely been ignored in previous studies of 
M13 because many (but not all) are found in the cluster core. The evidence of a small gap at 
V — I ^ 0.10 between the reddest BHB stars and the bright peak makes it worth considering 
whether they represent a separate population and what their evolutionary status is. For 
the purpose of this section, we restrict ourselves to comparisons with theoretical models, 
delaying additional discussion of the luminosity of the HB to §7.2.21 

In examining recent stellar models from different research groups, we found that there 
is still a great deal of variation in the morphology of evolution tracks for stars that start 
their HB lives with Tgfj ~ 10000 K. This is an importan t ques tion because it affects the 
interpretation of the reddest BHB stars in M13. ISweigartI (Il987| ) asserted that the relative 
importance of the helium-fusion luminosity {Lhe/L) is the main influence on whether the 
evolution is primarily blueward or redward. Stars that have strong hydrogen fusion shells 
significantly increase the helium core mass during the HB phase, leading to higher Lhs/ L and 
blueward evolution. As the helium becomes depleted, the importance of core fusion is reduced 
and stars tend to evolve redward, becoming more giant-like. For similar com p osit ions , some 
models predic t a modest blue ward loop before the approach to the AGB (e.g. ISweigartlllQSTl : 
Dormarull992l : lYi et al.lll997) for HB evolution start ing with Teg ~ 10000 K, while the most 



recent models (e.g. DSEP, 



Pietrinferni et a.1.1 120061 ) predict stars evolve redward (see Fig. 



fH\i . For the Teramo ( jPietrinferni et al.ll2006[ ) models, slightly cooler HB stars (ones starting 
among the reddest blue HB stars) have short blue loops. However, in the case of DSEP 
models, stars with —2 < [Fe/H] < all evolve strongly redward from the start, and remain 
close to the ZAHB during much of the HB phase. Consistent with Sweigart's discussion, 
core helium burning remains significantly stronger than the hydrogen fusion shell in DSEP 
models, even for stars starting quite close to the giant branch. If the DSEP models are good 
representations of actual evolution, even the evolved HB stars in M13 could be reasonable 
indicators of where the ZAHB is. Clearly there are significant differences between these 
models that may help us to identify important physics. 
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Although there are many small differences between the physics inputs in the DSEP and 
Teramo models, we believe the likely cause is diffusion. Diffusion was consciously disabled 
in the Teramo models, while the DSEP models stop diffusion in the outermost O.OSM© 
of stars and linearly ramp their diffusion algorithm to full strength at O.lOMp^ be l ow th e 



surface. As seen in more metal poor O.73M0 models with diffusion by lMichaud et al.l ( 120071 ). 
the introduction of diffusion produces redward evolution for a star that would otherwise 
evolve blueward. There is clear evidence of diffusive processes in the surface composition s 



of bluer HB stars in globular clusters (e.g. iMoehler et al.l l2003t iMoni Bidin et al.l |2009| ) 



but not among stars like the reddest BHB stars. Still, the surface layers are less likely 
to show chemical signatures due to near-surface convection. The DSEP models may be 
demonstrating that diffusion in the interior can be constrained by observations — processes 
such as rotationally-induced mixing could inhibit the action of diffusion. 

The numbers of stars can also constrain the models — because the early stages of 
core helium fusion take the longest, the initial direction of the evolution in color should 
have a big effect on how many stars are found to the red of peak PI. Fig. |22] shows a 
representative synthetic horizontal branch generated from DSEP web tool^ with a number 
of stars in the primary peak comparable to observations. We did not attempt to fully model 
the fainter parts of the primary peak because most of those stars are expected to exhaust 
their core helium before entering into the color range we are considering, becoming too bright 
to be considered HB stars (see §5.3.31) . Even so, the number of stars actually observed is 
considerably smaller than predicted by the synthetic HB simulation. The models predict that 
the red population (/ < 14.9) should be approximately 28% the size of the primary peak 
population, whereas the M13 population is only about 19%. Small variations in chemical 
composition (helium or heavy elements) also do not affect this conclusion. Whatever the 
reason, the DSEP evolutionary tracks do not seem to be accurately representing the evolution 
of stars from the primary peak. A small amount of blueward evolution would relieve the 
discrepancy. 

Teramo tracks produce a somewhat better fit to the M13 data, but the observations 
appear to need blueward evolution to continue at slightly higher temperatures to explain the 
appearance of the rather sharp red edge at {V — I) = 0.07 (pure redward evolution would 
tend to smear out such a feature in color). When the ZAHB and = 0.05 lines are fit to 
the magnitude extent of the primary peak, the reddest BHB stars and RR Lyrae stars fall 
near the Yc = 0.05 line, implying they are significantly evolved (see Fig. |23|1 . 



^http://stellar. dartmouth.edu/ models/shb.html We used "empirical" color transformations as these ap- 
pear to do a bet ter job than "synthetic" transformations in reproducing photometry in and around the 
subgiant branch (jPotter et aLlbOOTHSarajedini et al.1120071 ). 
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While the Teramo models come close to explaining the appearance of the sharp red 
edge of the HB distribution, the underlying question remains why there should be such a 
sharp edge for M13's population. Increased envelope mass, as well as increased helium and 
metal abundance are known to encourage blueward evolution. Helium enrichment has a 
significant effect on the HB luminosity and some effect on the evolution of red HB stars. 
For modestly helium-enriched {Y = 0.30) Teramo models in Fig. [2l] (consistent with the SY 
hypothesis), a large (and probably unrealistic) distance modulus of about 14.75 is required, 
and the evolution of stars at the red end of the bright peak PI goes too bright to explain 
the reddest BHB stars. In fact these stars can't be satisfactorily reproduced in the current 
helium-enhanced models. 



Pietrinferni et al.l ( 120091 ) recently computed HB models for "extreme" CNONa abun- 
dance mixtures in which the sum of CNO elements is approximately a factor of two higher 
at a given [Fe/H] than for a typical a-enhanced mixture. The mixture is intended to be 
representative of pollution resulting from intermediate mass AGB stars, and N is by far 
the most abundant heavy element due to nuclear processing. An extreme CNONa mixture 
could realistically produce blueward evolution for the dominant HB population. However, 
measurements of oxygen abundances among the cool HB stars in M l 3 all indica t e [0/ Fe] is 
normal for globular cluste rs and super-solar (jPeterson et al.lll995l ). ISmith et al.l ( 119961 ) and 
Cohen fc MelendezI (120051 ) also found that [(C-|-N-(-0)/Fe] did not seem to vary significantly 
within their samples of CN-strong and CN-weak giants, and the CNO elements do not show 
"extreme" enhancement (the average was abo ut 0.3 dex). Similar results have been found 
for relatively unevolved stars in other clu sters (ICarretta et al.ll20051) . although stars in NGC 



1851 show a significant 0.57 dex spread (lYong et al 



20091 



From the theoretical comparisons above, our preferred explanation is that the reddest 
BHB stars are stars with unenriched compositions that evolve somewhat to the blue after 
reaching the ZAHB. However, it must be admitted that the theoretical models disagree to a 
greater extent than we would like. We return to the discussion of these stars in §7.2.21 



5.3.2. Ultraviolet Bright Stars 



We next identified hotter ultraviolet-brig^ 
but optically id entifiable) from the study of 



it stars (th e brigh test cluster stars in U band, 
Zinn et al l Jl972l). ZNG 1 (Barn ard 29) is a 



post- AGB star JConlon et al.lll994l : iMoehler et"al1ll998[ [Thompson et al.l 120071 ) . and is the 
brightest star in the UV by almost 3 mag (see the UIT panel of Fig. [9]). 



In the lower left panel of Fig. El ZNG 2, ZNG 6, and ZNG 7 sit at the cool side of a 
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group of stars we identify as AGB manque stars, which are hot stars that are significantly 
brighter than t he ho t HB (see §5.3.3p . ZNG 2 (G 43) was observed spectroscopically by 
Moehler et al.l ( 120031 ). and its low gravity {logg ~ 3) clearly identifies it as being a post-HB 



star. ZNG 7 falls slightly bluer than the "knee" of the HB in the upper right panel of Fig. 
[TOl but is quite bright in the ultraviolet. ZNG 6 is the only one of the UV bright stars that 
has HST ACS observations, which again places it significantly brighter than the HB (see the 
lower left panel of Fig. [TOj) . 

In optical CMDs, ZNG3 and ZNG ^| are about a magnitude brighter than the knee of 
the HB, and part of a small group of stars seeming to parallel the HB. ZNG 3 and 4 are 
among the reddest objects with r eliable photonaetry d etected using UIT (see Fig. [9]). ZNG4 
was observed spectroscopically by lAmbika et al.l (120041 ). who identified chemical signatures of 
diffusion, which only shows up among blue HB stars with T^s > 11000 K. A likely explanation 
is that ZNG 4 {T^s = 8500 K) was once a hot blue HB star that evolved, and the chemical 
signature is leftover from the earlier phase. A significant convective envelope is not likely 
to appear until the star reaches lower Tes. We identify both stars as supra-HB stars, which 
probably trace evolution from midway on M13's HB, as discussed in the next subsection. 



5.3.3. Hot Post-HB Stars 



The origin and evolutionary behavior of the hottest HB stars has implications beyond 
stellar evolution — these stars contrib ute to, and perhaps dominate, the ultraviolet light 
from old stellar populations in galaxies (IDorman et al.lll993l ). Recently, iBrown et al.l (120081 ) 
used a STIS UV CMD of the dwarf elliptical M32 to study the main stellar contributors in 
the ultraviolet: the hot HB, the AGB manque, and the post-AGB. The authors found that 
the CMD could constrain the chemical evolution of the population, but that the post-HB 
evolution did not seem to be in good agreement with models. Although galactic populations 
provide better leverage on the shortest stages of post-HB evolution, a massive globular cluster 
with a simpler population and precise photometry should also constrain theoretical models. 

Using a combination of ultraviolet and optical CMDs, we attempted to identify all stars 
that have e volved away from the hotter parts of the H B. As discussed in §5.3. H theoretical 
HB models (IDotter et al.ll2007l : iPietrinferni et al.ll2006[ ) with no helium enrichment agree that 
stars originating from the primary peak in the HB distribution produce AGB stars exclusively 
— the evolution of the stars keeps them within about 0.2 mag of the HB until shortly before 
central He exhaustion. Core helium exhaustion occurs over a large portion of the star's 



For completeness, ZNG 5 is a non-member according to proper motions, and is not considered here. 
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core due to convection during helium fusion, causing a star to adjust its structure rapidly (in 
about a Kelvin-Helmholtz timescale) and move to the AGB clump. However, during the time 
leading up to core exhaustion, the core convection zone is decreasing in size, leaving behind a 
composition gradient. In typical AGB stars, almost half of the AGB lifetime is taken up while 
the new He fusion shell eats through this composition gradient. During this time the star's 
evolution pau ses, producing a fa irly well populated AGB clump about 1 V mag brighter 
than the HB (IFerraro et al.lll999l ). Figs. [T9] and [T8l show illustrative theoretical tracks and 
synthetic HB populations from DSEP and Teramo using canonical physics and reasonable 
choices for chemical composition (nearly primordial helium abundance, for example). The 
top rows in both plots show tracks for stars which evolve near the HB and have normal AGB 
phases including an AGB clump. 

The middle rows in Figs. Hn] and [TS] show stars that produce stars that spend part of 
their lives as somewhat bluer than average AGB stars, but do not have an AGB clump phase. 
For these stars, the evolution as the He fusion shell consumes the composition gradient take 
longer (more than a Kelvin-Helmholtz timescale) and occurs at higher surface temperature. 
Most of these stars can be put into one of two categories based on whether they have colors 
bluer or redder than the "knee" of the HB in optical CMDs {B — I ^ —0.1). One group of 
stars sitting approximately 1 mag above what would be the horizontal part of the normal 
HB (14 < B < 14.5 in Fig. [TOl V ~ 14 . 1) cor responds to "supra-HB" stars previously 



identified by IStrom et al.l ( Il970l ) and IZinnI ( 119741 ). This group includes ZNG 3 (AGB 75) 
and ZNG 4 (AGB 8), as well as AGB 27, 32, and 33 (detected in F160BW using HST, 
and the reddest post-HB star marked in the upper left panel of Fig. [H]), 50 (also detected 
in HST F160BW), 81, and two of the known BL Her pulsating variable stars (VI and 
V6; iKopacki et al.ll2003l ). Type II Cepheids (a group that includes BL Her stars) are only 
found among stell ar po pulations having a signifi cant blue HB com ponent, as was not ed by 
WallersteinI fll97oh and ISmith fc Wehlaul fll985h . As reviewed by IWallersteinI J^OoJ, the 
great majority of well-studied Type II Cepheids show no period change or increasing period, 
consistent with evolution toward the red and larger size. 

Even hotter post-HB stars are harder to identify in optical CMDs because of the steep- 
ness of the HB in magnitude. Relatively small color errors can cause normal HB stars to 
overlap stars that have evolved significantly in luminosity from hotter on the HB. When 
ultraviolet observations are used, these stars are much more easily separated by brightness 
(see Fig. |9]). So-called AGB manque stars fall in this category — they are stars that have 
evolved away from the HB that will not reach the traditional AGB because the star's en- 
velope does not have enough mass to produce a giant-like structure. Observationally, this 
definition is hard to use because it requires reliable knowledge of the star's future evolution 
in the CMD. However, theoretical calculations imply that there is a big change in track 



-27- 



morphology on the blue HB. The bottom rows of Figs. [T9] and [T8] show stars that would not 
have an identifiable AGB phase. 

The main difference in the models with total mass less than about 0.54 Mq at M13's 
metallicity is that hydrogen shell fusion is unable to provide much of the luminosity even 
during the core contraction following helium exhaustion. It appears that if Lh at maximum 
does not surpass Lue after the H fusion shell has reignited and stabilized, the star will remain 
at high temperature during the relatively slow helium shell adjustment. 

The supra-HB and AGB manque stars in M13 identified with crosses generally fall more 
than a magnitude above the faint envelope of HB stars. So for a fairly standard chemical 
composition, models predict that the only evolution tracks that produce hot post-HB stars 
come from bluer than some point between the two HB peaks in M13. Observationally, there 
is an apparent gap between two sets of "UV-bright" stars in the UIT CMD in Fig. O While 
individual optical CMDs in Fig. [TO] may leave the impression of color gaps between groups 
of hot post-HB stars, the union of the optical CMDs (see Fig. [23]) indicates that there is a 
thin, fairly uniformly populated band paralleling the HB from its bluest to its reddest colors 
with a couple of post-HB stars even closer to the AGB clump. There are a few stars in the 
ACS field without UV photometry that will partially fill the gap between the two groups of 
"UV bright" stars in the UIT CMD. As a result, we believe these stars are tracing out the 
slowest phase of post-HB evolution for stars that do not follow traditional AGB tracks. A 
relatively small change in envelope mass leads to a drastic change in the morphology of the 
tracks for these stars. Because M13's HB is well-populated near the end of the canonical HB 
(more so than many other clusters with extended blue HB tails), the cluster is providing us 
with a means of observationally "seeing" where the post-HB evolution pauses for a little. 

Independent of how believable the gap in the UIT CMD is, the usable portion of the 
UIT field does contain the majority of the hot post-HB stars, and does perhaps allow us 
to observationally identify the transition between an AGB manque track and an post- 
early AGB track. AGB manque evolution is expected to be largely in luminosity, with 
redward color evolution increasing for cooler, more massive HB stars. The red edge of the 
group (mi62o ~ y ~ ~^-'^) is thus a conservative upper limit for the HB stars producing 
AGB manquestars, and this limit is bluer than the blue end of the primary HB peak. The 
numbers of post-HB stars help to constrain the tracks further (see §6.21) . 

As can be seen from Fig. [2ll the shape of post-HB evolution tracks does not change 
drastically with a modest increase in helium abundance io Y = 0.30. Increasing helium 
abundance creates a larger luminosity gap between the blue end of the HB and the AGB 
manque stars, greatly increases the number of AGB manque stars r elative to HB stars , 



and produces bluer manque stars that are more concentrated in color ([Brown et al.ll2008l ). 
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Because models of M1 3's HB involving helium-e nriched stars predict that the bluest HB 
stars have Y ^ 0.38 (ID'Antona fc Caloil l2008bl ). these effects could be tested. Detailed 
simulations are beyond the scope of this study, but to first order, the number of hot post-HB 
stars is quite small compared to HB stars in the same color range (17 / 191 = 0.09 for 
^1620 — V < —1.9), which argues against large helium enrichment. 



6. Population Ratios 

One of the primary reasons for this study was to examine evolutionary timescales for 
stars in different life stages. We first recalculate the HB type {Nbhb — Nrhb) / {Nbhb + 
NvHB + Nrhb)- Here BHB, VHB, and RHB stand for stars bluer, within and redder than 
the RR Lyrae instability strip.) For M13, Nvhb = 9, and we have identified five possible 
RHB stars. Thus we have HB^ype = 0.976 ± 0.006. 



6.1. The Helium-Sensitive R Ratio 



The R = Nhb/Nrgb ratio compares the evol utionary lifetimes of HB and RGB stars 
and is a sensitive helium abundance indicator fe.g. ICassisi et al.ll2003r). Because of the claim 



that M13 stars may have high helium abundances ( Johnson fc Bolt^l998 : Caloi fc D'Antona 
20051). it is part i cularl y im portant to exa r nine t his population ratio. Following the definition 
of ISalaris et al.l (120041 ) and lZoccali et al.l (120001 ) . we used Vzahb (the magnitude of the zero- 
age HB) as the cutoff for the RGB sample. The determination of the faint limit for the 
RGB sample is important because the numbers of RGB stars rises quickly with increasing 
magnitude. 

To determine Vzahr, w e applied two methods. For the first, we followed the procedure 
of iRecio-Blanco et al.l (|2005[ ). using a template cluster of comparable metallicity and a well- 
studied RR Lyrae popula tion (used to define the reference HB level (Vr/j)) to determine a 
relative magnitude shift. iRecio-Blanco et al.l used NGC 1904 as the template for clusters 
near M13's metallicity after first determining (Vr/j) through a comparison with M3. In their 
comparison of photometry for NGC 1904 and M13 from the WFPC2 camera, they found 
that M13's sequences are 1.20 mag brighter in V than NGC 1904, with no relative shift in 
color. The corresponding {Vrr) value for M13 is therefore predicted to be 14.97±0.07, which 
results in Vzahb = 15.06. 



For the second method, we determined the ZAHB almost exclusively from data on the 
reddest blue HB stars, under the assumption that they are relatively unevolved and unen- 
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riched HB stars. (Based on additional arguments in §7.H we believe this is incorrect.) Our 
photometry does no t have adequate time -coverage for direct determination of (Vr/j). The 
value determined bvlKopacki et al.l ( 120031 ) (14.83 ±0.0 2) was se t by co mparison to photome- 



try by iRey et al.l (l200ll ) that was in turn calibrate d to iLandoltl (119921 ) standards. Zeropoint 
differences between our dataset and lKopacki et al.l should therefore be small. However, non- 
variable stars at the blue end of the instability strip have V = 14.90 on average. Using this 
as representative of the average HB level, Vzahb = 14.99. 

We find Nrgb = 483 using the first method, giving R = 795/483 = 1.65 ± 0.09 (error 
estimate from Poisson statistics). For the second method, Nrgb = 465 and R = 1.71 ±0.10. 
In both cases, M13's R value is about 3a higher than the theoretical calibration presented 
Salaris et al.l (120041 ) for Y = 0.245, and in agreement with their measurement oi R = 



m 



1.719 ± 0.197 using a considerably smaller sample. 

This does not account for the effects of HB morphology, however. Theoretical models 
uniformly predict that blue HB stars have longer evolutionary times than RR Lyrae variables 
or re dder HB stars , which would produce high R values in clusters with blue HB morpholo- 



gies. 



Zoccali et al.l ( l2000l ). for example, conducted an early examination of the effect of HB 



lifetimes on the R ratio, and ISalaris et al.l ( I2004J ) found that clusters with HBtype > 0.8 



show a larger spread around the mean than clusters with redder morphologies, which im- 
plied that the exact mo rphology of blue HB clusters might be influencing the R ratio. In 
Sandquist &: HessI (120081 ) . we described a method for correcting R for variations in HB life- 
times in a way that is based on observable quantities and is largely independent of chemical 
composition. The reader should see that article for more details, but we briefly summarize 
it below. 

According to theoretical models, HB lifetimes vary similarly as a function of stellar 
mass and effective temperature, although the absolute value of the lifetime does depend 
on composition (see the lefthand panels of Fig. [25]). Based on this behavior, we defined a 
weighting factor for each star in the HB sample: 



Wi 



tHB(JogTeff = 3.85) 



tHBi}ogT, 



eff) 



The largest portion of the variation in HB lifetimes occurs on the blue tail for Te// > 10^ 
K for stars with low-mass hydrogen envelopes. By choosing logT^ff = 3.85 (near the blue 
edge of the instability strip) as the normalization point, we can correct the blue HB lifetimes 
back to values representative of the more common variable and red HB stars, and the R 



(Salaris et al. 


2004; 


Zoccali et al. 


2000) 



'-eff 



IS 



shown in the righthand panels of Fig. [251 Although there is some variation in the weighting 
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factors with composition, they clearly describe the lifetime variation to first order, so that 
residual uncertainties are at the level of a few percent. Further, the weightings only remove 
the effects of position on the HB, and do not reference the absolute value of the HB lifetime, 
which depends on composi tion and physics in puts to the stellar evolution codes (e.g., the 
^^C(a;,7)^^0 reaction rate: ICassisi et al.l 120031 ). 



Although color- Tgj/ relationships remain imperfect, it is much more reliable to derive 
Tgff from photometry than it is to derive the stellar mass. For example, the Grundahl M-jump 



has b een found in many clusters with blue HB stars with T^ff > 11500 K (IGrundahl et al. 



19991 ). In addition, the canonical HB also seems to have a reasonably well-defined termination 
near 30000 K. For the pur poses of this study of M13, we have used the T ^jf determinations 
from iMoehler et al.l (120031 ) using Stromgren photometry. iMoehler et al.l also measured T^f f 
spectroscopically, and while the spectroscopic measurements are in good agreement with 
the photometric measurements, the spectroscopic measurements are subject to significant 
model uncertainties at the faint end and appear to deviate systematically from photometric 
determinations at the red end of the HB. These measurements for M13 do not go all the way 
to the faint end of the HB, so we have assigned T^ff = 31000 K to stars at the cutoff in the 



distribution (/ = 19.5) based on measurements of NGC 6752 (IMoehler et al.ll2000l ). which 
has similar metallicity and HB extent. 

As expected from the analysis above, the use of different sets of mo dels has a sma ll 
effect on the weighted HB star total. In the case of some model sets (DSEP, ISweigartlll987l ). 
we needed to extrapolate the weighting corrections in the range 4.4 < logTejj < 4.5 for an 
extremely small number of stars. Overall the r nodels agree that th e HB sample should be 
corrected downward by approximately 10%, with lCassisi et al.l (120041 ) (Y = 0.23, Z = 0.0006) 
and Sweigart {Y = 0.25, Z = 0.001) models giving A^^^ = 716, and DSEP (F = 0.248, 
[Fe/H] = —1.5) and Sweigart {¥ = 0.30, Z = 0.001) models producing a lower N'j^^ = 703. 
Thus, the corrected ratio is i? ~ 1.47±0.09 and 1.51±0.10 for our two methods of determining 
VzAHB- (The uncertainty introduced by the choice of models is significantly smaller than the 
Poisson uncertainties.) These values deviate from theoretical models for Y = 0.245 by about 
la. Based on these arguments, the R ratio app ears to rule out a global helium enrichment 
of AY = 0.04 at a 3a level (since dR/dY ^ 10; ICassisi et all 120031 ). 



The value of R produced by this analysis rests on the accuracy of the weightings used to 
correct the HB star total and on the determination of the faint cutoff for the RGB sample. 
If incorrect weightings were hiding a helium enrichment of Ay = 0.04 among M13 HB stars, 
they would need to have been overestimated by about 28%. We remind the reader that the 
weightings are corrections relative to HB stars near the instability strip and are independent 
of the absolute value of lifetimes, so we believe that an overestimate of this magnitude would 



- 31 - 



be difficult to produce. In addition, because we used weightings derived from models with 
0.23 <Y < 0.25, this would overestimate the weights for the bluest HB stars if those stars 
are actually helium-enriched, not underestimate them. If anything, our R measurement is 
biased slightly too high. 

To hide a helium enrichment through a systematic error in the RGB sample, it would 
be necessary to overestimate the number of stars, which would require a faint limit that was 
too faint. To produce a measured helium enrichment of Ay = 0.04, our faint limit would 
need to have been sXV = 14.79, or 0.20 mag brighter than our brightest estimate. Only a 
gross error in judging the level of the HB could produce an error this big, but the raajorit y 
of the RR Lyrae stars in the cluster dispute that, having ly) ^ 14.85 (IKopacki et al.ll2003l ). 
Even if the RR Lyrae stars are somewhat evolved, they are still fainter than the ZAHB 
level needed to produce a result of AY = 0.04. (The RR Lyrae stars will be discussed more 
below.) From other indications (such as AV^^), the HB of M13 is likely to be brighter than 
average for globular clusters if it is deviant at all. As a result, we don't see a clear reason to 
dispute the helium abundance implied by the R ratio. 

To avoid issues with systematic effects, it is wort h doing a relat i ve com pari son between 



M13 a nd M3. We assembled photometry for M3 from iFerraro et al.l (119971 ) and iRood et al. 



( 1999). along with avera ged photometry for the la rge population of RR Lyrae variables from 
Benko et al.l (120061 ) and ICorwin fc Carneyl (120011 ) . This sample completely covers a 7' x 7' 
area roughly centered on the cluster core. This contains 602 HB stars, which we subsequently 
corrected for lifetime effects. Because M3's HB morphology is much more horizontal than 
M13's, we have used the V — I color to determine lif etime corrections, in part due to questions 
about the calibration of the B data in this cluster (jValcarce fc Catelarul2008l ). For stars that 
had no I measurement, we determin ed an empirical tran sformation from B — V to V — I 
using nonvariable HB stars from the iFerraro et al.l ( 119971 1 data. For a handful of RR Lyrae 
variables that did not have average brightness information, we assigned the star a color equal 
to the average color of variables of the same type (RRab or RRc). The details are fairly 
unimportant because the normalization for the lifetimes is based on stars near the blue edge 
of the instability strip (near the middle of M3's HB) and because lifetimes are theoretically 
predicted to vary little except on the blue tail. For M3, the total lifetime correction amounted 
to just over 1%, with a corrected value of 594 stars. 

The cutoff magnitude for the RGB sample ((Vr_r)) for M3 can be determined quite 
accurately, so that we only need to worry about possible zeropoint difference s between the 
variable and non- variable star photometry. We used Vzahb = 15.70±0.03 from I Ferraro et al. 



(1997), which was also the source of the zeropoint for the HST observations of iRood et al. 
(Il999l ). The RGB sample brighter than Vzahb is 444 stars, which gives R = 1.34 ± 0.08. 
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According to this analysis, the M3 sample is consistent with having lower average helium 
abundance than M13, but the values differ by slightly more than the error bars and the 
implied difference in Y is only about 0.015. This is not enough to explain morphological 
differences in the CMD on the subgiant branch or HB. 

This is an interesting and somewhat surprising result that contradicts a number of pre- 
vious studies of M3 and M13 that indicated that M13 is enriched in helium compared to 
M3. For example, i n the ir discussion of globular clusters with multiple stellar populations. 



D'Antona fc Caloil ( l2008l ) identified M 13 as a cluster comp osed almost entirely of star s en 



riched in helium by about AY = 0.04. In a later simulation, iD'Antona fc Caloil (l2008b( ) also 
model M13's HB population exclusively with helium enriched stars, but in this case with 70% 
of stars having ranging from 0.27 < Y < 0.35 and 30% having Y = 0.38. For comparison, 
their interpretation of M3's populations involves 50% of the stars having a near-canonical 
value of y = 0.24 and 50% having 0.26 <Y < 0.28. Even if M13's R value is not corrected 
for lifetime variations, it still falls short of the value expected for the relatively small helium 
enrichment of AY = 0.04. We will discuss the cluster helium abundance in ^ 



6.2. The i?2 Population Ratio and Post-HB Evolution 

As the morphology of the HB becomes more blue, the evolutionary tracks of stars after 
the HB phase start to shift from traditional AGB tracks (starting in a clump at the faint 
end of the AGB and subsequently following a track paralleling the RGB) to abbreviated 
AGB tracks (having an evolutionary pause separate from the traditional clump, touching 
on the traditional AGB at moderate luminosities, and peeling away before reaching the tip) 
to manque tracks (retaining a surface temperature thousands of degrees higher than any 
part of the traditional AGB). M13 contains HB stars that clearly span this range according 
to the number of traditional AGB, supra-HB, and AGB manque stars we identified in §5. 
In view of the continuing difficulties in explaining the blue HBs of clusters like M13, we 
consider whether the AGB stars can reveal anything about the evolution or structure of the 
progenitor HB stars. 

At the simplest level, the R2 = Nagr/Nhr ratio compares the rela tive lifetimes in 



the AGB and HB evolutionary phases. IVargas Alvarez fc SandquistI (120071 ) found that for 



a sample of clusters with large samples of evolved stars (more than 200 HB st ars) that the 



popu lation ratio R2 dropped well below the theoretically predicted value of 0.12 ( ICassisi et al. 



20041 ) for clusters with a HB type > 0.8 (mostly bluer than the instability strip). When 
we include supra-HB and AGB manque stars in the AGB population, we derive a value 
i?2 = 90/795 = 0.113 ± 0.012. For the 4 clusters with the reddest morphologies (RBtype < 
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0.2), ISandquist fc Boltd f l2004h found {R2) = 0.106 ± 0.011. C l usters with 0.2 < HB 



type 



< 



0.8 had higher values (M5. 0.176 ± 0.018. ISandquist fc Boltd 120041: M55: 0.156 ± 0.023, 
Vargas Alvarez fc SandquistI 120071 ). Values for bluer clusters may be underestimated due 
to the difficulty of identifying AGB manque stars using only optical filters. Had we not 
included AGB manque and supra-HB stars in the AGB star counts, we would have calculated 
i?2 = 61/824 = 0.074 ± 0.009, completely consistent wi th the low values plotte d there for 
M30 and NGC 6752. In the blue HB cluster NGC 2808, ICastellani et~aD fl2006bh were able 
to identify manque stars, and found an R2 value similar to ours. With proper identification 
of the different types of stars, the R2 values are consistent with observational values for bluer 
clusters and with theoretical predictions. 

Based on the bimodal distribution of HB stars (approximately 47% of HB stars in the 
brighter peak, and about 28% in the fainter peak), we might expect to see evidence of an 
almost bimodal distribution of AGB stars with the manqueand supra-HB stars evolving from 
the bluest HB stars. The detected AGB manqueand supra-HB population is Nmanq/NACB = 
29/90 = 0.32 ±0.06 of the total population. If we assume as a first-order approximation that 
all AGB stars have equal lifetimes, then this implies that the manquestars originate in the 
bluest 32% of the HB stars (J > 17.75; V > 17.55). This excludes the bright peak and the 
majority of the intermediate HB population (see Fig. [T5|) . Taking into account that bluer 
stars have longer lives (and so are over-represented in the HB population), the red boundary 
for HB stars producing manquestars should be further to the red. When this is accounted 
for using the weighting factors discussed in §6.1^ the boundary falls at / ^ 17.3 {V ~ 17.15; 
^1620 — V ^ —3). Both imply the track morphology switches between the two maxima of 
the HB distribution, but closer to the fainter peak. 



7. Discussion 

M13 is distinguished by its horizontal branch, and its horizontal branch is notable for 
two reasons: the fact that the stars are almost exclusively bluer than the instability strip, 
and the fact that there is a secondary population of stars grouped near the extreme end of 
the horizontal branch. Below we discuss evidence bearing on each of these points and try 
to put M13 in context of other globular clusters. However, before we do, we would like to 
state that this is not intended as a complete survey of the subject, and we apologize for our 
limited ability to reference the many previous studies. 
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7.1. The Bimodal Horizontal Branch 



Independent of the position of the primary peak of the HB star distribution, the cause of 
blue HB tails has not been definitively identified. Some studies have identif ied relationships 



between HB tails and cluste r dynamical parameters su ch as central density (iBuonanno et al. 



19971 ) and total luminosity ( iRecio-Blanco et al.ll2006[ ). although neither of these can clearly 
explain differences between M3 and M13 because of their structural similarities (see Table 
1). While blue tails can plausibly result from a process (such as mass loss or chemical 
self-enrichment) producing a large dispersion in properties, M13 has a clear secondary peak 
in its HB distribution and this requires multiple star populations or the action of multiple 
physical processes. For example, theoretical models predict that varying amounts of mass 
loss on the RGB from star to star can produce a spectrum of outcomes from normal HB 
stars to early hot flashers (stars which leave the RGB before the helium flash, but ignite 
helium shortly afterward) to late hot flashers (stars which ignite helium on the white dwarf 
cooling curve) to helium white dwarfs (which never ignite helium). To produce a secondary 
HB peak in this way, some mechanism must be concentrating stars in the CMD, and it is 
difficult to see how this can be accomplished without almost complete loss of the hydrogen 
envelope for many stars on the RGB. On the other hand, the self-enrichment hypothesis 
explains peaks on the HB (not only at the blue end) via discrete populations of stars with 
different chemical compositions. A peak at the blue end of the HB can be accomplished 
with a population of stars having extreme helium abundances {Y > 0.35) because enriched 
stars leave the main sequence with lower total mass, and if the mass is low enough they will 
have almost no envelope by the time they reach the TRGB. He abundances for these hot 
stars probably can't answer the question directly (due to the action of diffusive and mixing 
processes), but correlated enrichment of h elium and carbon among the hottest stars in u Gen 
and NGG 2808 (IMoehler et al.ll2004j . 120071 ) is not predicted in the self enrichment picture. As 
yet, there is no clear way of distinguishing between these scenarios for EHB stars. However, 
clearer understanding of these stars may provide new clues, so we compare M13 with other 
clusters with blue HB tails below. 

Like some of the most massive clusters [NGG 2419, NGG 2808, uo Gen, and NGG 6273 
(M19)], M13 has a clear second peak at the blue end of the HB in optical filter bands. Even 
so, there is a range in the way the secondary peak appears. Some of t he best comparison s 



of different clusters (wit h HB s aligned in the op t ical) can be found in iPiotto et al 



f ll999h 



Dalessandro et al.l ( 120081 ) . and lSandquist fc Hesd (|2008[ ). Often there is a gap or edge feature 
in the HB star distributions as a function of an optical magnitude (My ~ 4.5, Mb ~ 4.2) 
and t his feature appear s to separate blue hook stars from the EHB. In M13 (and also NGG 



6752; ISabbi et al. 



2004J ). the secondary peak is almost entirely brighter than the position of 



the gap. The secondary peak in NGG 2419 reaches its maximum among the blue hook stars. 
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but appears to straddle the gap. The secondary peak in M19 straddles the position of the 
gap, but differential reddening prevents the clear identification of a feature, u Cen shows 
a sharp edge in the HB distribution at about the position of the gap in NGC 2419, with 
a much smaller fraction of stars just brighter than the gap. NGC 2808 on the other hand 
appears to have almost all of its hottest HB stars fainter than the gap, but its population is 
a considerably smaller fraction of its HB stars. Many other clusters with lon g blue HB tails 



do no t show a secondary maximum at the blue end of the HB, however. M54 ( iMomany et al. 
20041 ) and M15 appear to ha ve produced blue ho ok stars, but there is little or no sign of 



a secondary maximum. M80 (IFerraro et al.lll998l ) and M2 (IMomany et al.ll2004j ) have HBs 
with a similar extent to that of M13, but do not show maxima. So clusters are capable 
of producing EHB populations comprising from nearly up to 30% (for NGC 2419) of the 
total HB tally, with distinguishing features even among the very bluest HB stars in different 
clusters. 

As emphasized in models ( Brown et al.ll200ll ) and observations ( Momany et al. 2004), 
evolution of hot HB stars is mostly in luminosity and this can allow the evolved stars to 
masquerade as brighter but less evolved HB stars in optical photometry due to steep bolo- 
metric corrections. To better characterize the populations at the end of blue end of the 
HB, ultraviolet observations are needed. In near UV photometry (f/, F336W), the effects 
of luminosity evolution become better separated from color distribution. As can be seen in 
Fig. m far UV photometry (UIT B5, F160BW) mostly flattens the blue end of the HB in 
magnitude, and should allow the optimum separation of luminosity and color effects. Far UV 
observations are not very common, and it is even rarer to connect observations in different 
UV wavelength bands. 

A small number of clusters have archival data in the near and far UV as well as optical 
bands. HST WFPC2 data is available from proposals 5903 (M80; F. Ferraro PI), 6804 (NGC 
2808; F. Fusi Pecci PI), and 8709 (NGC 6752, M2; F. Ferraro PI). We processed images 
using HSTPhot in a manner similar to the M13 images. To compare different clusters, we 
identified stars at the red end of the Grundahl u jump in U band from HB morp hology 
or wi th the help of spectros copic information (NGC 6752, iMoni Bidin et al.l 120071 : NGC 
2808. iMoni Bidin et al.ll2009l ). The u jump appears to have a common T^s in all clusters 
(iGrundahl et al.lll999l ). and so we used it to register the CMDs by temperature. In add ition, 
blue HBs have a long segment at nearly constant U magnitude (IFerraro et al.l 119971 ) that 
includes the u jump, and in the F160BW band the u jump falls near the red end of a 
fairly flat segment reaching nearly to the end of the HB. Thus, this registration provides a 
convenient means of comparing the extent of blue HBs in color and magnitude. It also avoids 



the n eed to make large and uncertain corrections for reddening in U V bands (ICardel 



19891 ) or for the time dependence of WFPC2 throughput in the UV (iHoltzman et al. 



i et al. 



1995|). 
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Fig. |26] shows sets of CMDs shifted so that the red end of the u jump is at (0,0). 



A few conclusions can be drawn from simple comparisons. First, the large population of 
blue hook stars in NGC 2808 falls in a portion of the CMD {AU > 2, A{U -V) < -1) that 
is not occupied by more than a few stars in any other cluster. It is worth remembering that 
these are not the bluest stars in the A F160BW, A(F160BW - V)) CMD. They are fainter 
(A F160BW> 0.5) and slightly redder than the blue end. This is the origin of the "blue 
hook" monik er, and it is though t to be the result of a helium-rich atmosphere resulting from 
flash mixing (IBrown et al.ll200ll ). 



M13, M2, and M80 have similar morphologies at the blue end of the A F160BW, 
A(F160BW — V)) CMD, though the distributions of stars differ. The HB sequence dips 
faintward in the range —2.4 < A(F160BW — l^)) < —2 before rising again up to (F160BW — 
V)) = —3. Most of the stars in the rising portion are also in the Momany U jump (see Fig. 
|27|) . The identification between the upturn in F160BW and the Momany U jump is clearest 
in M13 because of the large population of stars at the end of the HB. Some stars are found 
to the red in the {U, U — V) CMD for M2 and M80, although these may be blends with main 
sequence stars in V bancl^. 

The correlation between stars in the Momany jump and the upturn in F160BW is 
not perfect, but it suggests that there is a coherent group of stars to be found on the 
extreme blue HB. The group does not fall at the precise end of the HB though, as we 
find in each cluster small numbers of stars that are bluer and fainter in the F160BW and 



U CMDs. This interpretation differs from that of iMomany et al.l (120041 ). who stated that 
the HB sequence appeared to make an almost discontinuous jump in color (~ 0.3 mag) 
at constant U magnitude, and extending faintward in U (by ~ 0.5 — 0.7 mag) at nearly 
constant color. The color histograms make it clear that the three clusters have very different 
distributions, with M2's distribution declining toward the blue end, M80 with a broad but 
evenly spread group, and M13 with a well-defined peak and a large fraction of the stars in 
the Momany jump. 



While no spectroscopic data exists in the literature for M13 stars in this range. iMoehler et al 



( 120031 ) examined cooler HB stars (8000 K < T^s < 21000 K) and found evidence of strong 
helium depletion and iron enrichment for Tes ^ 1 2000 K. We can get ad ditional guidance 
by contrasti ng spectroscopic resu lts for NGC 6752 (iMoni Bidin et al.ll2007l) against those for 
NGC 2808 dMoehler et al]l2004l ) and u Cen dMoehler et al.ll2007[ l. NGC 6752, u Cen, and 



^Main sequence stars contribute very little light in UV bands, and so blends shift almost horizontally to 
the red in the CMD. Because the F160BW— F color has a longer wavelength baseline than U — V, blends 
don't shift stars as much. 
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NGC 2808 all have noticeable vertical jump features in [U, U — V) CMDs, but the most ob- 
vious jumps occur at larger T^s in NGC 2808 and u Cen. In NGC 6752, there is consistently 
helium depletion in the atmospheres of stars hotter than the Grundahl jump at T^s ~ 11000 
K, including the hottest stars observed. In NGC 2808 and u Cen, stars with near-solar and 
super-solar helium abundances (and concurrent carbon enrichment) appear at temperatures 
above those seen in NGC 6752 (Tcs ~ 31000 K). Both spectroscopic signatures are consis- 
tent with the predictions of "late hot flashers" : stars that ignite helium on the white dwarf 
cooling curve, and initiating convective mixing of the envelope. 



We are led to the same conclusion reached by lMomany et al.l (120041 ): that the stars in the 
second U jump are likely to be "early hot flashers". In this case, the existence of a working 
hydrogen-fusion shell is thought to inhibit convective mixing into the outer envelope during 
core helium ignition. As a result, no clear spectroscopic signature is expected and none has 
been found to date. Because heliu m enrichment of the envelope plays a role in producing 



blue hook stars ( Brown et al.lboOll ) by reducing atmospheric opacity shortward of the Lyman 



limit (and reducing redistribution of flux to longer wavelengths), we exar nined the data for 



spect roscopically studied stars in the Momany jump within NGC 6752 (IMoni Bidin et al. 



20071 ). We found no correlation between measured helium abundances [which covered the 
range —3.26 < \og{nHe/nH) < —1.58] and U magnitude, which seems to confirm that the 
envelope is still dominated by hydrogen. The scatter in helium abundance appears to be 



interesting, however — when lMoni Bidin et al.l ( 120091 ) combined their measurements for M80 



and NGC 5986 stars with those from NGC 6752, they found consistent helium depletion 
for stars hotter than the Grundahl jump at Tgfj ~ 11000 K, but also found stars with 
larger depletions and generally larger star-to-star scatter in two temperature ranges (13000 
K < Teff < 18000 K and 25000 K < Teg < 31000 K). 

Both of the helium depletion features appear to be associated with the U jumps, al- 
though the details probably differ. For the Grundahl jump, stars that are hotter than the 
jump are uniformly brighter in the Stromgren u filter than a canonical ZAHB. By con- 
trast the stars associated with the second U jump in M13 are brighter than the majority 
of HB stars of similar (but higher or lower) temperature, as seen in the U and F160BW 
filters. Observationally, the radiative levitation of heavy elements and the downward dif- 
fusion of helium almost certainly must be involved because all stars with Tefj < 11000 K 
(excepting blue hook stars) are found to have some degree of surface helium depletion, but 
the mechanics are still being debated. The transition at 11000 K i s associated with the 



near-disappearance of the surface convection zone, but recent models (IMichaud et al.l 120081 ) 
indicate a need for a low-mass mixed layer near the surface to moderate the build-up of 
metals. There may be reason to attend to the effects of rotation because a discontinuity in 
average rotation speeds also appears at around this temperature (with hotter stars largely 
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being slow rotators; iPeterson et al.lll995t lBehrll2003[ ) 



Based on the evidence above, we considered whether progenitors of extreme HE stars 
such as early and late hot flashers can be identified on the RGB. Early hot flashers are 
theoretically expected to only leave the RGB near the tip, while late hot flashers could leave 
at lower luminosities. While M13 appears able to produce both kinds of hot flashers, the 
second U jump stars (our early hot flasher candidates) are much more abundant. Thus, we 
would expect that clusters with similar HB morphologies (like NGC 6752, M2, and M80) to 
show symptoms of this near the TRGB if anywhere. 



Strangely, ISneden et al.l (120041 ) found that M13's known "super 0-poor" stars ([0/Fe] 
< —0.4) have high Na abundances and most appear very close to the RGB tip (My < —2.3), 
indicating the possib le exposure of heavily processed gas at the surface. For perspective, 
Carretta et al.l (120061 ) summarize literature data for the well-kriown a nticorrelation between 
[0/Fe] and [Na/Fe] for stars in 20 clusters, and ICarretta et al.l (120081 ) present new spectro- 
scopic data for 19 clusters. The extent of the 0-Na anticorrelation varies from cluster to 
cluster, but it probably requires nuclear processing under conditions that cannot be produced 
in gas that can be mixed to the surface in low-mass stars. This has led to the supposition 
that it results from pollution by previous generations of more massive stars. Even so, super 
0-poor stars are quite rare and, to the best of our knowle dge, have not been s een among 
relatively unevolved turnoff and subgiant branch stars (see ICarretta et al.l 120061 ) . although 
this could be due to difficulties in detecting low oxygen abundances among faint stars. 

In addition, the cluster NGC 2808 has a large population of blue hook stars and 
seemingly few or no early hot flasher candidates, and contains super 0-poor stars as well 
(ICarretta et al.ll2006l ) but these are fairly uniformly spread down the RGB (although those 
authors did not observe stars at the tip of the RGB), u Gen also has a large population of 
blue hook stars, and appears to have a population of super 0-poor stars (INorris fc Da Gosta 
19951 ). In the cases of NGG 2808 and u Gen though, there is clear evidence of both self- 
enrichment and multiple populations, which may complicate the interpretation. 



Examination of bright giants in NGG 6752 (lYong et al.ll2003l : iGarretta et al.ll2008) have 



revealed no super 0-poor stars. In addition, Mg isotope ratios in M13 (ISneden et al 



20041 ) 



and NGG 6752 are correlated with the O depletions among bright giants, but the conditions 
necessary to process Mg are not expected in RGB stars based on current nuclear reaction 
rates. The sample of iGarretta et al.l (120081 ) hints at the possibility of smaller populations of 
super 0-poor stars in clusters like NGG 3201 and M5, neither of which have extreme HB 
stars. 



Though there is strong evidence of mass loss effects on the extreme BHB, there is still 
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little to connect this directly to the properties of stars still on the RGB. One outstanding 
question that deserves theoretical attention involves the position (redder than the extreme 
end of the HB) and clumping of the early hot flasher candidates in the CMD. Even if helium 
is enriched overall in the bluest HB stars, some non-standard physics is probably required to 
produce their observed characteristics because helium-enriched EHB stars should be fainter 
than ones with primordial helium abundance. We can suggest two types of observations that 
might help to clarify the situation. First, in clusters with stars concentrated near the blue end 
of the HB, stars at the TRGB may have chemical peculiarities similar to those in M13 if they 
are preparing to have early hot helium flashes. If found, this could result from excessive mass 
loss or from highly helium-enriched stars, but it would connect RGB and extreme HB stars. 
Second, in clusters with predominantly red HBs, spectra for blue HB stars in a limited range 



of Tp ff near the instability strip may reveal their initial helium abundances ( jVillanova et al. 



20091 ) . This would be a direct observational test of whether helium enrichment has produced 



any blue HB stars. While there is circumstantial evidence of helium enrichment in clusters 
like NGC 6388 and NGC 6441 (IBusso et al.l 120071 ). spectroscopic proof should be possible. 



7.2. "Second Parameter" Effects on the HB 

Because the R ratio seems to rule out previously suggestions of helium enrichment in 
M13, we discuss the literature on M13 to see if the observations can be reconciled. The main 
evidence pointing toward a helium enrichment falls into several categories: 1) morphology 
of the HB, 2) luminosity of the HB and the RR Lyrae variables, and 3) morphology of the 
subgiant branch. 



7.2.1. HB Morphology 

M13 and M3 have long been a "second parameter pair" based on the distribution of 
stars on their HBs — the small difference in [Fe/H] (the "first parameter") between the two 
clusters is unable to explain the much bluer HB stars in M13. Before discussing further, 
we emphasize that we are mainly considering the shift in the CMD position of the "most 
representative" HB stars in the two clusters, and not the differences in the shape of the HB 
distributions. Our goal is not to model the HB morphology in detail. M13 has an unexplained 
bimodal HB distribution as well as a significant number of stars between the two peaks, but 
arguably the brighter peak (15 < < 15.5; slightly bluer than the instability strip) is the 
most representative based on star numbers. By contrast, M3 has a u nimodal distribution 



with the peak of the distribution falling within the instability strip (jValcarce &: Catelan 
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20081). 

Earlier studies ( iFusi Pecci et al.lll993l : iBuonanno et al.lll997l ) have used {B — V)peak (the 
dereddened color of the peak of the HB star distribution) as a way of describing the most 
representative HB stars in clusters, although this indicator starts to lose sensitivity on the 
steep blue tail. We used our earlier results from histogramming the I magnitudes of HB 
stars, and then translated the resulting position to the {B — V) co lor using a fiduci a l line . 
We find Ipeak ~ 15.25 and {B — V)peak ~ —0.05. For comparison, iBuonanno et al.l (Il997| ) 
give {B — V)peak = 0.30 for M3. If the clusters have nearly identical chemical compositions, 
this corresponds to a difference of about O.OGMq in HB star mass (Teramo models). Unfor- 
tunately, this does not transform to a similar mass difference at the cluster turnoff if common 
red giant mass loss formulas are applied because they are nonlinear with star mass. If M13 
is helium enriched compared to M3 by AY = 0.05, the mass difference is virtually the same. 
Helium enrichment in M13 does not remove the need for a difference in mean mass between 
its HB stars and those in M3. 

The mass difference can be explained with reasonable assumptions in both the At and 
Ay hypotheses, and we refer the reader to §1 for a brief summary of recent attempts to model 
M13's HB. In the At hypothesis, M13's stars are older than M3's. In the AY hypothesis, 
helium enrichment allows less massive stars to leave the main sequence earlier, but a large 
amounts of mass loss is still necessary for each star. While this may be correct, the picture is 
not fully motivated by physical reasoning. As discussed in §5.3.1^ the sharp red edge of the 
bright peak and the population of redward-evolving HB stars may provide new constraints 
on fits in the AY hypothesis. We encourage synthetic HB studies that pay closer attention 
to these features. 



7.2.2. The Luminosity of the HB 

The luminosity of the HB is most strongly affected by metallicity and helium abundance 
{6VHB/Sy ~ —4 mag), and is insensitive to age. However, without accurate and independent 
measurements of distance and extinction, the luminosity must be determined relative to other 
cluster stars. Most commonly, the reference landmarks are found among MS stars (such as 
the turnoff or a fainter point identified with the help of the turnoff) or RGB stars (like the 
RGB bump). These methods have generally indicated that M13's HB is luminous compared 
to other clusters, but the reference points usually have the disadvantage of having their own 
dependencies on composition and age that complicate the int erpretation. In the discussion 
below , we will focus on the AY = 0.04 Jjohnson fc Boltelll998l ) and At = 1.7 Gyr feev et al. 



200ll ) hypotheses for explaining differences between M3 and M13, but we must consider the 
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evolution hypothesis (see §1) at the same time. 



When the cluster turnoff is used as a reference, i t introduces a dependence on age 
{SVTol^t ~ 0.075 mag / Gv r: ICaloi &: D'Antonal 120051) and increases the dependence on 
helium {SYtoI^Y ~ 1-5 mag; ICaloi fc D'Antonall2005h . We determined Vto = 18.59 ± 0.05, 
which gives AVr^^^ = 3.69 ± 0.07 if we apply our value for Vhb a t the blue edge of the 



instability strip. This value is large r than any prev i ously quoted (see lJohnson fc Boltelll998 



for a summary), and if we use the iKopacki et al.l ( l2003l ) value for (Vrr), the value would 
be higher still (3.76 mag). The larger values of AV^^^ depend entirely on whether the RR 
Lyra es and "reddest blue HB " sta rs near the i nstab ility strip (see §5.3. ip are evolved or not. 
Both lJohnson &: Boltd ( 119981 ) and iRey et al.l (l200l[ ) observed much smaller samples in M13, 
and an examination of their CMDs shows that together they only observed two "reddest 
BHB" stars among the 55 we identified — un derstandable cons idering that most of these 
stars are in the core of the cluster. As a result, IJohnson fc Boltd quot e a value for Vhr thai 
is 0.19 ma g fainter than ours. We also note that like the studies of iJohnson fc Boltd and 
Rey et al.l . our photometry for the HB and turnoff regions are on a consistent zeropoint. We 
see no reason to dispute earlier values of AV^^^ for M3, given that M3's HB is well-populated 
on either side of the instability strip. So if M13's reddest HB stars are not evolved, it has 
an extreme value for AV^ff (almost 0.2 mag larger than that of M3), and this is consistent 
with both the At and AY hypotheses. To find the cause (the HB, TO, or both), we need 
other reference points. 



Ferraro et al.l (Il997| ) conducted a similar comparison between M3 and M13 using HST 



U observations. In U, the subgiant branch and a section of the blue HB become flat, making 
them excellent magnitude references. The authors found no significant difference between 
the AU^Q^ derived for the two clusters. While this does not directly bear on the luminosity 
of the HB near the instability strip, it is an indication that the bluer HB of M13 has a normal 
luminosity compared to M3. 

The red giant bump can also be used as a referenc e point, dependent on heli um abun- 



1.5 



mag; 



Riello et al. 



2003 



dance {SVbump 
{SVbump/St ^ 0.035 mag / Gyr; ICaloi fc D'Antona 



Caloi fc D'Antonal |2005|) and age 



20051). The observed value for the 



difference AV^^^ ca n be put in context using the catalogs of iFerraro et al.l (1999) and 
Di Cecco et al.l (120101 ). To ensure a reliab le comparison, we calculated the correction from 
Vhb to VzAHB = 14.96 according to the iFerraro et al.l prescription. We therefore derive 

among those for clusters of similar metallicity (for ex- 
ample, M3 has —0.23 ± 0.07). However, their tabulated value for Vz ahr is 0.14 mag fainter 
than ours — this is partly a reflection that the photometry they used (jPaltrinieri et al.lll998[ ) 
lacked stars near the instability strip, and partly that they used synthetic HB calculations 
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to derive the level of the ZAHB relative to the HB stars they did observe. Using the fainter 
ZAHB, M13's value for AV^^^ becomes significantly different from M3's value, and falls at 
the low end of the distribution for similar clusters. So if the RR Lyraes are showing us the 
true ZAHB level, neither age nor helium changes are necessary but can't be ruled out because 
the expected changes are la level. If the ZAHB has a fainter level, rather large differences 



are needed {AY = 0.06 - 0.09, or At = 3.5 Gyr). The tabulation of iDi Cecco et all (120101 ) 
measured values for both M3 and M13 using homogeneous photometry and a template-fitting 
method for determining Vzahb for troublesome clusters like M13 that have very blue HBs. 
Their AV^^^ values are in agreement to within the measurement errors!^ 



The brightness of the AGB clump is fairly independent of metallicity ( iCastellani et al. 



19911 ). but when the magnitude difference with the HB is formed, the de pendences on 



helium abundance and RGB progenitor mass almost completely disappear (jPulond Il992 



Cassisi et al.l 120011 ). Although M13 has a sparsely populated AGB clump, there is a fairly 
clear grouping of 12 stars at = 14.19. Using the bright ZAHB level from the previ- 
ous paragraph, we have AV^p^ = —0.7 7. This is slightly smaller than the value for M3 



(-0.88 ± 0.10) from iFerraro et al.l (11999f ). but only about Icr different. If the fainter ZAHB 
level is applied, the M3 and M13 values agree to within about 0.03 mag. The results might 
also be affected by likely lower mean HB mass in M13 compared to M3, which affects how 
the clump is populated and would tend to reduce the luminosity of the clump. We consider 
this weak evidence in favor of the fainter ZAHB level. 

The HB luminosity can be constrained using the end of canonical HB as a reference in 
clusters with long blue tails. Like NGC 6752, M13's HB appears to terminate essentially at 
the end of the canonical HB, and we identified blue hook candidates (in §5.2p and second 
U jump stars (in §7.ip that mark points in the HB near the end. Models show that the 
difference in magnitude between the termination of the blue HB and the blue end of the 
instability strip grows substantially larger with increasing helium abundance, mostly because 
the horizontal part of the HB gets brighter. The faint end of the HB is fairly insensitive to 



ic 



Caloi fc D'Antonal (|2005| ) also used the position of the RGB bump relative to the turnoff (AV^™^) as a 



potential helium indicator, with larger magnitude differences implying higher helium. In their comparison, 
they found that M13 had a magnitude difference that was 0.14 ± 0.09 larger than that of M3, which they 
admitted had relatively low significance. Because our photometry of the bump and main sequence constitutes 
a homogeneous dataset with a much large r number of stars, we checked their quoted bump and turnoff 
magnitudes (derived from the ph otometry of Paltrinieri et al.lll998r) . We find there is a rather large difference 
in color (0.08 mag) between the iPaltrinieri et al.l photometry and ours, but the V magnitudes are almost 
identical to ours. On their own, the RGB bum ps of the M13-M3 pair do not give a strong constraint on 
helium abundance, although lCaloi fc D'Antonal showed that clusters with predominantly blue HBs had larger 
Ay^Q™^ than clusters that have strong populations in the instability strip. 
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helium abundance because the hydrogen-rich envelopes of the stars are almost gone. If the 
end of theoretical ZAHBs are shifted in magnitude to fit the end of M13's HB, the bright end 
is consistent with primordial helium values (see Fig. [25]) . While the fit is understandably 
uncertain due to the gap between the HB stars and the blue hook stars in optical filters, 
the Teramo evolutionary tracks do go nearly parallel to the faint envelope of HB stars (see 
Fig. |23|) . The fit seems to rule against helium enrichment among stars within and redder 
than the primary peak. Evolutionary effects among the reddest HB stars would tend to bias 
toward an indication of helium enrichment, so this result seems fairly robust. Spectroscopic 
measurements of stars near the red end of the HB in M13 and NGC 6752 also support the 
idea that the outer envelopes have not been significantly pro cessed through the C NO cycle. 



and (in NGC 6752) still have primordial helium abundance (IVillanova et al.ll2009[ ). 



Alternately, if the ZAHBs are fit to the red end of the HB, helium-enriched ZAHBs are 
far too red at the faint end, lending additional credence to models with canonical composition. 
However, current models of the bluest HB stars do not fully incorporate important physics 
such as diffusion and radiative levitation, and these effects seem to be responsible for the u 
jump and the intra-peak HB stars that are brighter than canonical models. As a result, fits 
to the HB should be treated somewhat skeptically. 

Finally, we consider the distance moduli implied by the HB (specifically the RR Lyrae 
stars) and by the TRGB. Both of these features are in common use as standard candles, and 
a comparison of the resulting values might also give us clues on whether the HB is unusually 
bright. The absolute magnitude of th e TRGB in / is a l most independent of age and chemical 



composition. From the calibration of iBellazzini et al.l ( l2004l ). we calculate Mj^^^ = —4.07 



which gives {m — M)i = 14.35. The distance modulus could be smaller than this if by chance 
M13 does not contain RGB stars very near the flash stage. To evaluate this possibility, we 
can use a binomial distribution to calculate the probability that at least one star is within 
a certain magnitude range of the TRGB. In our cluster-wide sample, there were 271 stars 
brighter than the RGB bump (J < 13.6). Using theoretical RGB luminosity functions that fit 
M13's population (see §1]), stars brighter than that level have an approximately 0.4% chance 
of being within 0.01 mag of the TRGB, so that a calculation using the binomial distribution 
predicts a 70% chance of having at least one of the 271 stars within that interval. The 
probability rises to about 90% that at least one star is within 0.05 mag of the TRGB, so 
this is probably a minor source of error. 

All of the brightest RGB stars are known to be semi-regular variables including the 
brightest (Vll) , having an arnplitud e of about 0.13 mag in V and a timescale for variability 



of about 30 d (IKopacki et al. 



20031 ). As a result, the variability of t he brightest gia nts is 



a significant source of error in our determination of Itrgb- Though iKopacki et al.l made 
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observations over 70 d (more than two cycles for most giant variables), it is possible they 
did not sample the full amplitude of variation. Our CFHT V me asurement of V ll (which 
was used above) was 0.15 mag brighter than the average given by iKopacki et al.l . while our 
KPNO observation was only about 0.06 mag brighter (barring zeropoint differences). If we 
can assume that the / observations are off by similar amounts, this gives (m — M)/ ^ 14.5. 
The range of variation for other bright RGB stars (V17, V24, V38, and V42) overlap that of 
Vll, and all but one give corrected distance moduli within 0.03 mag of the Vll value. V17 
gives a significantly brighter value (14 .37), but it also has the largest amplitude (0.38 mag) 
of the red giant variables tabulated by lKopacki et al.l . We believe the agreement among 4 of 
the 5 variables is sufficient to identify the TRGB, which in turn implies (m — M)o = 14.48 
with a very conservative uncertainty of 0.10 mag. 



Although there is still considerable discus sion of the RR Lyrae My — [Fe/H] relation, 
recent consensus views (e.g. ICatelan et al.ll2004j) implv a me an RR Lyrae value M^^ = 0.62 
for M13's metallicity. With (Vrr) from the iKopacki et aD feo03h observations of M13 RR 
Lyraes, this gives (m — M)v = 14.21 . The disagreement with the TRGB distance modulus 
can again be relieved if the RR Lyraes in M13 were significantly evolved or were helium 
enriched (Ay ^ 0.05). If we use the TRGB distance modulus and My^ to calculate the 
brightness of the HB, we find Vhb = 15.07. This is fainter than all of the "reddest BHB" 
stars but approximately the same brightness as the blue end of the primary peak. 

Once a gain, a compa r ison with M3 is useful, thanks to that cluster's huge RR Lyrae 
population. lJurcsik et al.l ( 120031 ) identified 4 groups of RR Lyraes with different brightness 
and hght curve parameters, with the most abundant group of 50 stars having (V^) = 15.67. 
Their most evolved group (32 stars), which the authors hypothesized were evolving redward, 
had {V) = 15.53. The difference between (V) for the most evolved RR Lyraes in M3 and that 
of M13's RR Lyraes is most consistent with magnitude differenc es derived from comparing 
points like the RGB bump and AGB clump (IFerraro et al.l 119991 ). 



Evolution is also expected to affect the pulsational properties of stars within the instabil- 
ity strip because the mean density of evolved stars is lower than for unevolved stars, resulting 
in larger periods. We could expect to discern evolutionary effects in several ways: based on 
observations derived solely from the light curves (such as the Bailey period-amplitude di- 
agram), period-luminosity relationships (with luminous RR Lyraes having longer periods), 
and color distributions (with evolved stars likely to be spread more evenly through the 
instability strip due to their short evolutionary timescale). 

If variability information (such as period and variation amplitude) can give unambiguous 
evolutionary information, it should be preferred because such measurements can generally 
be made to higher precision than can determinations of average magnitude or color. Fig 
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shows V amplitude and period data for M13 RR Lyrae stars from iKopacki et al.l (|2003[ ) in 
the Bailey diagram. The single known RRab clust er variable (V8) has a n even longer period 
than the "well-evolved" stars in M3 identified by ( ICacciari et al.ll2005l ). Because V8 is also 
the reddest RR Lyrae and evolutionary tracks in this part of the CMD are expected to slope 
brightward toward the AGB, it would not be surprising if it showed an evolution signature. 
In M13, the three BL Her stars (Population II Cepheid variables) fall approximately in the 
part of the CMD where stars that have evolved from the blue HB tail are predicted to pass, 
and their periods (all greater than a day) also support evolved status. However, V8's {V) is 
fainter than most of the RRc variables. Given the amount of scatter seen in larger samples 
of RR Lyraes (see below), this can't be considered evidence. 

RRc variables are more abundant in M13, so an evolutionary signature among these 
stars would be more compelling. Four of the RRc variables have periods consistent with the 
"well-evolved" RRc stars in M3, but with somewhat lower amplitude (one exception falls on 
the M3 relation). The remaining 4 (V7, V31, V35, V36) have periods consistent with the 
regular variables in M3 but again with low amplitudes. There is no evidence of separation in 
average V magnitude between the two groups however, weakening the idea that evolutionary 
effects can be identified solely from the periods of the RRc variables. 

To put M13 in better context, we assembled literature data for other Oosterhoff group 
II clusters (hereafter, OoII). For a few in this group, the instability strip is heavily pop- 
ulated, leading to an expectation (based on evolution timescales) that a large proportion 
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Silbermann fc Smithlll995l ). These clusters have relatively low BBtype values (< 0.8), 
and NGC 2419 and M68 have significant red HB populations as well. Both of these factors 
imply that there are likely to be unevolved stars in the instability strip. This group of OoII 
clusters, however, is outn umbered by clusters like M13 that have a lmost exclusively blue HB 
stars. For example, M2 (ILazaro et al.l l2006l : iLee fc Carneyl Il999l ) is a massive cluster that 
has more RR Lyrae stars than M13, but still has a very small fraction of its stars in the 
instability strip. For comparison, we compiled data for clusters with only a few RR Lyraes, 
under the hypothesis that these stars are more likely to have evolved from the blue HB0. 
The Bailey diagrams {V amplitude versus period) for the different samples are shown in Fig. 



"References for RR Lyra es in "BHB Clusters" are: M9 Jciement fc Shelt w] llQQOl) 
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Cohen fc Matthewslll992l): NGC 28 8 kaluznv et al.l Il997h : NGC 5897 iwehlaul ll99ol : Iciement fc Rowel 



20011) : NGC 5986 (lAlves et al.ll200ll) . 
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M13 RRc variables fall in parts of the Bailey diagram that are commonly occupied by 
variables in other OoII clusters. In the BHB cluster group and most of the other clusters, 
a significant fraction of the RRc variables fall near {Ay, log P) = (0.45, —0.41). This group 
appears largely independent of cluster metallicity. The main exceptions are NGC 5286, 
which is more metal rich than M3 and has an analogous group that is much closer in period 
to M3 RRc variables, and M2, which doesn't have a clear concentration of RRc variables 
in the Bailey diagram. The second identifiable group of M13 RRc variables overlaps with 
groups in the BHB Cluster sample and in co Cen, and has logP ~ —0.5 and Ay < 0.5. From 
these diagrams, we conclude that there are not clear evolutionary signatures that be inferred 
from the Bailey diagram alone. 

Strangely, there are proportionately few short period, low amplitude RRc variables in 
the "protoypical" OoII clusters M15 and M68. In addition, both of these clusters have 
populations of double-mode RRd varia bles that largely reside in a small rang e of color 
between the RRab and RRc variables (jSilbermann fc SmithI Il995l : IWalken 119941 ) . Such a 
distribution is unlikely if most of the stars are evolved BHB stars — evolution accelerates 
as a star moves redward toward the AGB. Conversely, we find that the shortest period 
(logP < —0.46) RRc stars are the bluest RR Lyrae stars in OoII clusters. Many of the 
known non-radial pulsators in OoII clusters (iKopacki et al.ll2003l : lOlech et al.lll999l : iKopacki 
200ll ) are also found within this group. 



At best, we find that Bailey diagrams for OoII RRc stars provide some information 
about average colors. This is somewhat useful for examining the effects of evolution within 
individual clusters because some color distributions (especially ones biased toward the red) 
are incompatible with model predictions. However, the RRc don't appear to follow patterns 
implied by M3, in which "well evolved" stars appear shifted to greater period. 

In summary, the majority of the evidence from photometric indicators implies that 
the M13 RR Lyraes (and nonvariable stars of similar color) are significantly brighter than 
the ZAHB level. The most notable exception involves the comparison to the RGB bump, 
while the pulsation properties of the RR Lyraes themselves are somewhat ambiguous. If the 
reddest HB stars are discounted as significantly evolved, then M13's blue HB appears to fall 
at approximately the same luminosity level as M3's and the evidence for helium enrichment 
or age differences from this is removed. 
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1.2.3. Subgiant Branch Morphology 



As mentioned earlier, one of the more unusual aspects of M13's CMD is the steeper 
slope of its subgiant bran ch in comparison to M3 's. This was inte rpreted as result ing either 
from helium enrichment (jjohnson fc Boltdll998l ) or greater age (iRey et al.ll200l[ ). Both of 
the cited papers used the so-called "horizontal method" to compare the ages of clusters with 
other clusters or isochrones. Higher helium abundance or age tends to shorten the length of 
the subgiant branch, so that if the turnoffs are aligned in color and fainter main sequence 
points (usually 0.05 mag redder than the turnoff ) are subsequently aligned in magnitude, the 
relative position of the giant branch reveals the difference. It has been verified repeatedly 
that there is a difference between M3 and M13, but the cause is unclear. 

We can attempt to look at the problem from a different angle by realizing that age and 
helium enrichment affect the absolute colors of the turnoff and red giant branch in different 
ways. An increase in helium abundance reduces the opacity of the stellar envelope, making 
surface temperature higher at both the turnoff and on the red giant branch. On the other 
hand, increased age makes both the turnoff and red giant branch redder, although the effect 
on the giant branch is very small. 

Because reddening and metal abundance differences have large effects, reliable compar- 
isons of absolute colors can only be done wh en these are w ell-determined — M3 and M13 
are probably the best example of such a pair. IStetsonI (Il998l ) discusses them in this respect, 
but we expand on the arguments here. The reddenings for the two clusters appear to be 
small and very sini ilar (M13 probably with the larger reddening by AE{B — V) < 0.01; 
Schlegel et al.lll998[ ) and their metallicities ap pear to agree to within 0.1 dex (with M13 be- 
ing the more metal poor; ISneden et al.ll2004j ). One additional benefit of using the M3/M13 
pair is that the effects of the reddening and metallicity differences on the turnoff should 
partly cancel, with the metallicity differences leading to an expectation that M13's turnoff 
should be bluer than M3's by about 0.01 in B — I, and less in other optical colors. 

A comparison in absolute colors also requires datasets for which the calibration can be 
done uniformly for both clusters. While we are unable to do this because our deep exposures 
of M3 and M13 were not taken on the same p hotometric night, i t has been done in the most 
recent studies, although using different colors. IJohnson &: Boltd (Il998l ) used V — I and found 
that M13 was bluer at the turnoff and giant branch by similar amounts (0 .01 — 0.02). This 
implies a tiny helium enrichment [AY ^ 0.01) at most. On the other hand. iRey et al.l (120011 ) 
used the B—V color and found that the turnoff colors were virtually identical, although M13's 
re d giant branch was still significantly bluer. We examined the standard star photometry 
of [stetson] (120001 ) for the two clusters (8 November 2007 update), and found that in the 
V — I color Ml 3 is bluer than M3 on the giant branch and at the turnoff (in agreement 
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with iJohnson fc BoltdllQQSl ). while in B — V or B — I colors M13's turnoff is significantly 
redder (even accounting for small differences in reddening) and the giant branches have nearly 
identical colors (see Fig. l5U]l l^. Based only on the Stetson standard stars in B — V and B — I 
colors, an age difference appears tenable. With the V — I color included, neither helium nor 
age differences seem capable of explaining the observations because neither is expected to 
affect one color differently than others. We conclude that current photometric datasets do 
not paint a consistent picture, and may still be influenced by subtle systematic effects. To 
make this a stringent test, effort should go into deep and carefully calibrated photometry 
using filters with a wide wavelength baseline. The B — I color still appears to be a good 
choice. 



8. Conclusions 

To our minds, some of the most important questions regarding M13 remain in dispute. 
One question that we have reopened here is whether the reddest of the blue HE stars in 
M13 are significantly evolved (and whether they are therefore good representatives of the 
brightness of the horizontal branch). The weight of the observational and theoretical evidence 
leans toward the idea that they are significantly evolved, and that the red edge of the primary 
HB population is a decent indicator of the true HE level. 

The distribution of stars on the horizontal branch in M13 is complex, and the most 
notable questions regard 1) how the color of the primary peak in M13 could have been shifted 
so far relative to M3's when the gross composition of clusters appear nearly identical, and 2) 
how a large fraction of M13's RGE stars become blue stars near the end of the canonical HE. 
Our examination of the luminosity function shows little sign that a large fraction of stars 
leave the bright RGE before h aving a core fiash, in agreer nent with the massiye clus ter NGC 



2419 f lSandquist fc Hessll2008h but not with NGC 2808 fISandquist fc Martelll2007f ). We do 
not find any clear evidence of helium enrichment among the stars of the dominant (redder) 
HE population, and in fact, the helium abundance indicator R and the relative brightness of 
the HE argue against significant enrichment. The HE and RGE stars (and different subsets of 
these) do not show significant signs of radial segregation within the cluster. The small color 



""^^While the Stetson standard stars in M3 and M13 are rigorously tied to the system of Landolt standard 
stars, the observations for the two clusters were taken under varied conditions. M3 and M13 can be observed 
on the same night using t he same equipin ent, and so it should be possible to get good relative photometry 
for the pair. We note that Stetson ( 19981) reported that M13 had a bluer turnoff than M3 by 0.014 mag in 
B — I (about 0.04 mag when reddening was accounted for) for images of the two clusters taken on the same 
night under photometric conditions using the same equipment. 
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difference between the main sequence turnoff and the giant branch of M13 (in comparison 
to M3) remains unexplained, but careful examination of the absolute colors of both clusters 
would provide a new test. 

Our thorough search of M13's HB population has revealed second U jump and blue 
hook stars that imply that many of these stars have very low-mass hydrogen-rich envelopes. 
Far UV observations show that many of the stars in the second U jump are more luminous 
than stars with similar colors at the end of the HB. The reason is u nclear, however. Spectro- 



scopic data on similar stars in NGC 6752 (IMoni Bidin et al.l 120071 ) indicate that the excess 



brightness is not related to enhanced atmospheric helium abundance, so further study is 
required. 

Spectroscopic measurements may help to clarify our understanding of extreme HB stars 
in a number of ways. We particularly encourage studies of: stars n ear the red end of th e 



HB in MIS where helium abundances can be accurately determined (jVillanova et al.ll2009l ): 
stars in the extreme HB to look for signs of unusual Mg abundances (a species that appears 
to be minimally affected by diffusion) that could connect them to giant stars; O, Na, and Mg 
for stars at the red giant tip of other clusters to determine whether they are super 0-poor; 
relatively unevolved turnoff and subgiant stars in M13 and NGC 2808 to search for large O 
depletions ([0/Fe]< —0.4) and check whether this is the result of external pollution or not. 
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100 -100 



X (arcsec) 

Fig. 1.— Outlines of the observed HST fields for M13. The HST proposal 5903 field with a 
dotted border, the proposal 8278 field with a solid border, and the proposal 8174 fields with 

a dashed border. (Three overlapping fields were observed as part of proposal 8174.) The 
ACS fields of the two WFC chips are shown with a dark sohd line (WFCl on the left). 
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Fig. 2. — Photometric re siduals between our calibrated photometry and the standard mag- 
nitudes of Jstetsonl[2nnnl ) in the sense of this study minus Stetson's. Stars measured on chip 
11 of the CFHT field are shown with open circles, and stars from chip 12 are shown with 
crosses. 



-62 - 




Fig. 3. — Photometric residuals between our calibrated CFHT photometry and the ACS 
photometry cahbrated to the standard system (in the sense of CFHT minus ACS) . 



-63- 



0.2 
-0.2 



Ix ' ' ' ' 1 ' ' ' 'x 1 ' ' 'x' 1 ' ' M 

— ^ Xo< — 

- X ^ ^ X X 

XX*^ X !( X X X X xX 


1 1 1 1 1 1 'x 'x 1 X x' 1 1 

X X X X 

" X xX % X "'xx XX Xx X - 
- Xx %|x X x^ 
Xx X x## x„ \,xX X 


- XX ''S.x xXX X ^ x "x " 

^ 1 , , , , 1 , , , , 1 , , . , ' , X , ' " 


r*<x " x^S^**" • X X : 
X xX^xX'x xX 

- x ,« f ^ 

x xxx xxXx X 

- Xx ^ 

" ' , , X , 1 X , , , . , , , - 



15 16 17 18 19 14 16 

V I 



-0.2 



1 J Ijj 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

— >*< X — 
X X xX^^ ^ X - 

- x/x€*«x' x'- 


1 1 jx 1 1 1 ^1 1 1 1 1 1 1 1 1 1 1 1 
X X X - 
_ X X 

x»<*'^>5x'*«*' X- 
X X „x fefx<xx X x^ 


^'^"^iMtP^ ^ 

: X :Skx*x'< >. ; 

1 ^ ^ /x^^ .xj 
, 1 X, , , , 1 , , , , 1 , , , , 1 


Xx xx^ Xx X - 

X X X X !< X X _ 
X 

— XX — 

, 1 , , ^ , ^< , , , , 1 , , , , r 



-0.5 0.5 1 -0.5 0.5 1 

V-I V-I 



-0.2 





M ' ' ' 1 ' ' ' ' 1 

''Xx xx" XX X* 

X ^ xX X X 


1 1 1 1 1 

« - 




-^x -xW^'xx--^ 

X x^!>:x*^ 

X >« X 

X? X 

1 . . - 1 


— ^ 

- 

X 

, , , , r 



-0.5 0.5 

V-I 



Fig. 4. — Photometric residuals between our calibrated CFHT photometry and the WFPC2 
photometry (F555W and F814W filters) calibrated to the standard system (in the sense of 
CFHT minus WFPC2). 
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Fig. 5. — Phot ometric r esidua ls between our calibrated KPNO photometry and the standard 
magnitudes of IStetsonI (120001 ) . 
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Fig. 6. — Photometric residuals between our calibrated KPNO and CFHT photometry (chip 
11) in the sense of KPNO minus CFHT. The sample has been restricted to stars brighter 
than the main sequence turnoff (with the exception of faint HB stars). 
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Fig. 7. — Photometric residuals between our calibrated KPNO and CFHT photometry (chip 
12) in the sense of KPNO minus CFHT. The sample has been restricted to stars brighter 
than the main sequence turnoff (with the exception of faint HB stars). 
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Fig. 8. — Color-magnitude diagrams used in the identification of bright RGB stars. RGB 
stars have AGB stars have A, and known field stars have x symbols. 
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Fig. 9. — Ultraviolet color-magnitude diagrams used in the identification of HB stars. Vari- 
able HB stars have * symbols, AGB manque and supra-HB stars are + symbols, "UV-bright" 
stars have Q symbols, and the hottest HB stars have y symbols. In the UIT panel, the 
V magnitude comes from the images with the highest available resolution (ACS, CFHT, or 
KPNO). 
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Fig. 10. — Optical color-magnitude diagrams used in the identification of HB stars. Symbols 
are the same as in the right panels of Fig. [HI with the exceptions that AGB stars are shown 
with A symbols, and unevolved non-variable HB stars are shown with □ symbols. 
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Fig. 11. — Comparison of the cumulative luminosity function of M13 RG B stars (solid line) 



with theoretical predictions from a Victoria- Regina model [dotted line; IVandenBerg et al. 



20061 ) with [Fe/H]= —1.41 and [a/Fe] = +0.3. The theoretical model has been shifted 



horizontally to match at the tip of the RGB, and vertically normalized to the RGB sample 
at / = 13.7, just brighter than the RGB bump 
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Fig. 12. — Comparison of the combined cumulative luminosity function of M13 and M5 RGB 
stars (solid line) with t heoretical predictions from the Victoria- Regina model (dotted line; 
VandenBerg et al.l 120061 ) with [Fe/H]= —1.41 and [a/Fe] = +0.3. The theoretical model has 
been shifted horizontally to match at the tip of the RGB, and vertically normalized to the 
RGB sample at / = 13.7, just brighter than the RGB bump 
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Fig. 13. — Combined color-magnitude diagram for bright giant stars (AGB stars are identi- 
fied with crosses) from all photometry sources. 
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Fig. 14. — The differential luminosity function of M13 RGB stars with the bump and possible 
hump identified. Error bars are based on Poisson statistics. 
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Fig. 15. — /-band Distribution of horizontal branch stars in M13. Top panel: the total 
sample. Middle panel: stars within 1 half-light radius (94"). Bottom panel: stars outside 1 
half-light radius. 
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Fig. 16. — Cumulative radial distribution of star populations. Top panel: red giant stars 
{solid line), and horizontal branch stars {dotted line). Bottom panel: horizontal branch stars 
divided by magnitude into faint peak (/ < 18; dotted line), intermediate stars (16.25 < I < 
18; solid line), red HB stars (/ < 14.9; dashed line, and combined RGB/HB population {dark 
solid line). 
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Fig. 17. — Zoomed color-magnitude diagram from ACS data. Symbols are as in Fig. [TOl 
The ZAHB [solid line) and the Yc = 0.10 locuses [dotted line) are shown for a-enhanced 
Teramo HB models having Z = 0.001 and Y = 0.246, with [m - M)qoq = 14.55 and 
E(F606iy - F8UW) = 0.02 assumed. The AGE clump falls just above the bright end of 
the plot. 
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Fig. 18. — Comparison of grou nd-based CFHT and KPNO photometry (r > 200") with 
Teramo (jPietrinferni et al.ll2006l ) HB models. In all panels theoretical values are shifted by 
14.45 in V and 0.02 in 5 — / to fit the envelope of HB stars. Left panels: Comparison with 
evolution tracks having [Fe/H]= -1.62 and [a/Fe] ~ +0.4 ([M/H] = -1.27, Z = 0.001). 
The ZAHB and the Yc = 0.10 locuses are shown with dotted lines. The tracks are for 0.4912, 
0.50, and 0.52 Mq {bottom panel), 0.54 and 0.56 Mq {middle panel), and 0.58, 0.6, and 0.63 
Mq {top panel). 
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Fig. 19. — Comparison of ground-based CFHT and KPNO photometry (r > 200") with 
DSEP HB models. In all panels theoretical values are shifted by 14.28 in V and 0.02 in 
5 — / to fit the envelope of HB stars. Left panels: Comparison with evolution tracks having 
[Fe/H]= —1.5 and [a/Fe] = +0.4. The ZAHB and the = 0.05 locuses are shown with 
dotted lines. The tracks are for 0.50 and 0.52 M© {bottom panel), 0.54 and 0.56 M© {middle 
panel), and 0.58, 0.6, and 0.65 M© {top panel). Right panels: Example synthetic horizontal 
branch populations having the same initial composition and a mass dispersion of 0.01 M©, 
The number of synthetic HB stars roughly corresponds to the number in the total sample 
(see Fig. [T5l) . Also included {solid line) is an isochrone for the same composition and an age 
of 14 Gyr. Crosses show our mode fits to MS data, and individual stars that represent the 
mean HB line. 
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Fig. 20. — (/, V — I) color-magnitude diagram for HB stars from all photometry sources. 
RR Lyrae stars (shown with * symbols) are found within the ACS field, and the photometry 
used is not averaged over the oscillation period. 
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Fig. 21. — HR dia Krams for HB stars from theoretical models. Top pa nel: Models from 
Potter et all ^20071 ) {soM line) for Z = 0.00107, [a/Fe] = +0.4 and from IPietrinferni et al. 



(j2006|) {dashed lme) ior Z = 0.001 ( 0.52 0.54, 0.56, 0.58, 0.6, 0.65 M© for both sets.) Bottom 
panel: Models from lPorman et all Jl993h f or [Fe/H] = -1 48, [0/Fe] = 0.63 {soM Ime; 0.52, 
0.54, 0.56, 0.59, 0.61, 0.65 Mq) and from IyI et aD f |l997h for Z = 0.001 {dotted Une; 0.52, 
0.54, 0.56, 0.58, 0.60, 0.64 Mq). 
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Fig. 22. — Comparison of observations {upper panels) with DSEP models {lower panels; 
assuming (m — M)/ = 14.28) of the primary peak on M13's HB. Left panels: {I, V — I) 
color-magnitude diagrams. Right panels: /-band distributions. 
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Fig. 23. — Comparison of mastered photometry with Teramo (jPietrinferni et al.l 120061 ) HB 
models. In all panels theoretical values are shifted by 14.55 in V and 0.01 in \^ — / to fit 
the envelope of HB stars, and have [Fe/H]= -1.62, Y = 0.246, and [a/Fe] ~ +0.4 ([M/H] 
= —1.27, Z = 0.001). The ZAHB and the Yc = 0.05 locuses are shown with dotted lines. 
Evolution tracks are shown for 0.57, 0.60, and 0.62 Mq {upper left panel), 0.54 and 0.56 Mq 
{lower left panel), and 0.4912, 0.50, and 0.52 Mq {right panel). 
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Fig. 24. — Comparison of mastered photometry with Teramo (jPietrinferni et aLll2006l ) hehum 
emiched [Y = 0.30, Z = 0.000924) HB models. In all panels theoretical values are shifted 
by 14.75 in V and 0.01 in V — / to fit the envelope of HB stars. The ZAHB locus is shown 
with a dotted hne. Evolution tracks are shown for 0.57 and 0.61 Mq {upper left panel), 0.54 
and 0.56 Mq {lower left panel), and 0.481, 0.50, and 0.52 Mq {right panel). 
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Fig. 25. — HB lifetimes as a function of T^ff from DSEP f lDotter et al.l 120071 ) models for 
(from top to bottom) [Fe/H] = —1.0,-1.5, and —2.0 ([a/Fe] = 0.2), and Sweigart models 
Y = 0.30,0.25, and 0.20 {Z = 0.001). Left panels: lifetimes in Myr. Right panels: lifetimes 
normalized to the stars reaching the zero-age HB at logTg// = 3.85. (For the DSEP model 
with [Fe/H] = —2.0, the most massive models do not reach that temperature, and were 
normalized at logT^fj = 3.89.) 
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Fig. 26.— CMDs of 5 globular clusters with blue HB tails in the F160BW and F336W (cal- 
ibrated to U) filters from HST. Crosses show stars identified at the red end of the Grundahl 
u jump, and all CMDs have been shifted so that the bluest of these stars is at (0,0). 
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Fig. 27. — Similar to Fig. |26l but zoomed on the blue end of the HB. A show are stars 
that were identified in the jump feature in the F160BW filter, x were ones identified in the 
F336W ([/) filter, and Q o^i^s that were identified in both. 
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Fig. 28. - — Comparison of grou nd-based CFHT and KPNO photometry (r > 200") with 
Teramo (jPietrinferni et al.ll2006h ZAHB models for Y = 0.246 {solid line), 0.30 {dotted line), 
and 0.35 {dashed line). All ZAHB curves have been shifted by 0.02 in 5 — J, but by 14.45, 
14.40, and 14.30 mag respectively in V. Blue hook candidates are shown with y symbols. 
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Fig. 29. — C omparison of period- amplitude values for RR Lyrae stars in M3 ( iCacciari et al. 



20051 ). M13 (IKopacki et al.ll2003l ). and other Oosterhoff II group clusters. For M3, stars 



identified as "evolved" (long period / high amplitude) are shown with x, and small amplitude 
variables are shown with A (among the RRab stars, these are suspected Blazhko variables). 
For other clusters, fundamental mode pulsators are shown with •, first harmonic pulsators 
have as double-mode pulsators have -|- (and are plotted with the dominant first overtone 
period), and non-radial pulsators have *. Solid lines are mean relations f or regular M3 
v aria bles, and dotted lines are mean relations for "evolved M3 variables" ( ICacciari et al. 



20051 ). In (jj Cen, stars are segregated by magnitude: for (y) > 14.6, RRab stars have Q 
and RRc stars have •; for 14.4 < (y) < 14.6, RRab stars have • and RRc stars have Qs 
for ly) < 14.4, all RRab stars have x. 
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Fig. 30. — Stetson standard star values for M13 and MS at the cluster turnoff in different 
optical colors. The vertical hues are the median color values of stars within about 0.2 mag 
of the cluster turnoff in brightness. 



Table 1. Characteristics of Ml 3 and M3 







M13 


M3 


[Fe/H]cG97 ^ 




-1.39 


-1.34 


Ra (kpc) ^ 




21.5 ±4.7 


13.4 ±0.8 


Rp (kpc) 




5.0 ±0.5 


5.5 ±0.8 






0.62 ±0.06 


0.42 ±0.07 


^ (deg)*" 




54 ±5 


55 ±2 


(kpc km s~^ 




-376 ± 145 


705 ± 123 


Etot (10^ km2 s- 




-0.476 ± 0.088 


-0.649 ± 0.022 


MvY 




-8.66 


-8.65 


logpo (/^VoPC"^' 




3.33 


3.51 






2.11 


1.56 


(pc)'= 




1.7 


1.3 


rc (") 




46 


26 


r/, (pc)'' 




3.5 


3.7 


TH ar 




94 


73 


d (kpc)'^ 




7.7 


10.4 


E(5 - vy 




0.017 


0.013 



Carretta fc GrattonI (119971 ) 



TOinescu et al.l fll999h 



McLaughlin &: van der Marell (120051 ): power-law 



models were used for Mvf values 



Schlegel et all (119981 ) 



Table 2: Description of HST data 



Proposal ID Principal Investigators Filters 



5903 Ferraro F160BW,F255W,F336W,F439W,F555W 

8174 van Altena F555W,F785LP 

8278 Bailyn F555W,F814W 

10775 Sarajedini F606W, F814W 



Table 3. Star Cross-Identifications 



ID 


L ID 


K ID 


B ID 


J. ACtj^ 


Other Notes 






V5 


806a 














RRl; 


blend in KPNO 


V7 


344 


432 











RRl 




V8 


206 


389 











RRO; 


reddest RRLyr 


V9 


806b 














RRl 




V25 

















RRl; 


evolved? bad sampling?; SIMBAD L630 ID is wrong (nearby RGB) 


V31 


807 





1043 








RRl; 


av affected by blending; BARN201 (blend) 


V34 


918 














RRl 




V35 


571 














RRl; 


gap in ACS 


V36 

















RR2 





Note. — Table |3] is published in its entirety in the electr onic edition of the Ast ronorn i cal Jo urn al. A portion is shown here for 



guidance regarding its form and content. Sources for IDs: L: iLudendorffl (Il905l ): K: iKadlal ( Il966l ): A: lArpI (Il955f l: B: iBarnardI (Il931l ) 



^UIT Photometry Flag. 0: No UIT detection; 1: Optimal photometry with no likely blending effects; 2: Photometry probably 
minimally affected by neighbors; 3: Photometry contamination of uncertain magnitude; 4: Photometric contamination certain; 5: 
Blending of two UIT sources. 



Table 4. HB Star Photometry 



ID 


Aa 


AS 


B 


V 


/ 


Source toi620 


(71620 Flag 


F160BW 




F255W 


0'255 


F336W 




PpM 


V5 


57.711 


1.414 


15.252 


14.690 


14.308 


ACS 









16.180 


0.042 


15.487 


0.004 




V7 


-52.404 


-67.743 


15.283 


14.944 


14.726 


ACS 









16.196 


0.181 


15.352 


0.005 


87 


V8 


-106.089 


25.345 


99.999 


14.922 


14.428 


ACS 

















99 


V9 


58.994 


3.471 


15.003 


14.745 


14.363 


ACS 





18.010 


0.118 


16.231 


0.071 


15.297 


0.003 




V25 


14.847 


-5.506 


14.753 


14.589 


14.368 


ACS 





17.240 


0.100 


15.590 


0.042 


14.954 


0.003 




V31 


57.728 


78.754 


99.999 


14.890 


14.565 


ACS 

















99 


V35 


91.867 


-28.731 


15.070 


14.833 


14.579 


ACS 









16.487 


0.191 


15.143 


0.003 




V36 


-0.226 


10.760 


14.834 


14.676 


14.349 


CFHT 





17.209 


0.083 


15.999 


0.060 


15.257 


0.004 





Note. — Table mis published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form 
and content. S 



Table 5. AGB Star Photometry 



ID 


Aa 


AS 


B 


V 


/ 


Source 


TO1620 


O'l620 


Flag F160BW crieo F255W ctsss 


F336W (7336 PpM 


1 


-612.315 


557.213 


17.108 


17.398 


17.725 


KPNO 


12.83 


0.03 


1 




2 


-403.576 


148.486 


15.225 


15.395 


15.524 


CFHT 


12.84 


0.03 


1 


99 


3 


-233.816 


112.551 


16.444 


16.659 


16.834 


CFHT 


13.42 


0.04 


1 




4 


-222.956 


29.858 


14.766 


14.144 


13.336 


CFHT 









99 


5 


-184.247 


134.084 


15.793 


16.004 


16.194 


CFHT 


12.96 


0.02 


1 





Note. — Table [S] is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form 
and content. 



Table 6. RGB Star Photometry 



ID 


Aa 


AS 


B 


V 


/ 


Source F255W 0-255 F336W 0336 


PpM 


1 


-707.261 


7.295 


14.373 


13.365 


12.236 


KPNO 


99 


2 


-556.761 


489.164 


15.380 


14.515 


13.533 


KPNO 


99 


3 


-543.549 


-318.220 


15.149 


14.304 


13.285 


CFHT 


99 


4 


-468.740 


-80.657 


15.734 


14.968 


14.024 


CFHT 


99 


5 


-439.069 


44.185 


15.543 


14.739 


13.785 


CFHT 


99 



Note. — Table |6] is published in its entirety in the electronic edition of the Astronomical Journal. A portion 
is shown here for guidance regarding its form and content. 



Table 7. Bright Nonmember Star Photometry 



ID 


Aa 


AS 


B 


V 


/ 


Source 


PpM 


K178 


-785.649 


-294.677 


14.410 


13.740 




PM 





K180 


-763.792 


-269.343 


12.280 


11.070 




PM 





K185 


-725.361 


-382.757 


15.380 


14.700 




PM 





K190 


-702.285 


186.036 


15.784 


14.926 


14.078 


KPNO 





K196 


-664.843 


556.797 


15.468 


14.764 


14.019 


KPNO 


12 



Note. — Table [7] is published in its entirety in the electronic edition of 
the Astronomical Journal. A port ion is shown here f or guidance regarding its 
form and content. I D No tes: L: iLudendorfi (Il905 ): K: IXadlal |l966h : CM: 
Ciidworth & MonetJ |l979h : SA, SB: ISavedofj (|l956h : F: this paper. 
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Table 8. HB Star Distribution 



Sample 


/ < 16.25 


16.25 < / < 18 


/ > 18 






/pi 


// 


/P2 


/Pl — /P2 


Total 


365 


197 


222 






0.47 ±0.02 


0.25 ± 0.02 


0.28 ±0.02 


0.18 ±0.02 


r < rj2 


25 


10 


21 






0.45 ±0.07 


0.18 ±0.05 


0.38 ±0.06 


0.07 ±0.09 


r < Tc 


82 


36 


61 






0.46 ± 0.04 


0.20 ±0.03 


0.34 ± 0.04 


0.12 ±0.05 


rc/'2. <r <rh 


145 


80 


99 






0.45 ± 0.03 


0.25 ± 0.02 


0.30 ±0.03 


0.14 ±0.04 


rh <r < 3.5rh 


157 


94 


87 






0.46 ± 0.03 


0.28 ± 0.02 


0.26 ±0.02 


0.21 ± 0.04 


r > 3.5r^ 


38 


13 


15 






0.58 ± 0.06 


0.20 ±0.05 


0.23 ±0.05 


0.35 ± 0.08 
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Table 9. Kolmogorov-Smirnov Test Results for Comparisons with Victoria-Regina RGB 

Models 



M13 Sample Combined M13 + M5 Sample 

I Cutoff D P N I- Itrgb Cutoff D P N 



13.4 


0.080 


0.10 


232 


3.3 


0.036 


0.46 


545 


13.2 


0.078 


0.19 


193 


3.1 


0.040 


0.44 


460 


13.0 


0.068 


0.45 


159 


2.9 


0.052 


0.23 


397 


12.8 


0.098 


0.16 


130 


2.7 


0.039 


0.69 


329 


12.6 


0.110 


0.13 


109 


2.5 


0.048 


0.56 


273 


12.4 


0.102 


0.25 


97 


2.3 


0.046 


0.69 


234 


12.2 


0.093 


0.43 


86 


2.1 


0.049 


0.72 


200 


12.0 


0.126 


0.19 


73 


1.9 


0.073 


0.30 


171 


11.8 


0.118 


0.40 


55 


1.7 


0.097 


0.12 


148 


11.6 


0.139 


0.34 


44 


1.5 


0.083 


0.40 


114 


11.4 


0.130 


0.66 


30 


1.3 


0.087 


0.48 


91 


11.2 


0.161 


0.58 


22 


1.1 


0.091 


0.59 


70 



